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I.  INFORMATION  SYSTEMS 


A.  Project  6151.  PATTERN  RECOGNITION 

Principal  Investigator:  T.  M.  Cover 

Staff:  R.  King,  A.  El  Gamal,  W.  Rogers,  K.  Jarett, 

J.  van  Campenhout , P.  Algoet 

1 . Objective 

The  purpose  of  this  project  is  to  investigate  problems  of  in- 
formation processing  and  pattern  recognition. 

2 . Current  Status  of  Work 

a.  Measurement  Selection  (T.  Cover,  J.  van  Campenhout) 

Virtually  any  ordering  of  the  probabilities  of  error  P(S) 

as  a function  of  the  2°  subsets  S cr  (X  ,X  , . . . ,X  } is  possible.  As 

■”  l ^ n 

a result,  it  can  be  demonstrated  that  no  nonexhaustive  measurement-se- 
lection procedure  exists  for  finding  the  best  k-element  subset  of  mea- 
surements. This  work  is  being  extended  to  prove  that  arbitrary  P(S) 
assignments  are  achievable. 

b.  Kolmogorov  Complexity  and  Shannon  Entropy  (T.  Cover, 

C.  Leung) 

The  Shannon  and  Kolmogorov  notions  of  complexity  can  be 
shown  to  be  essentially  the  same  by  placing  the  algorithmic  problem  of 
Kolmogorov  in  a probabilistic  framework. 

c.  Estimating  the  Entropy  of  English  (T.  Cover,  R.  King) 

A certain  gambling  estimate  of  entropy  is  asymptotically 
consistent.  Experiments  reveal  that  the  entropy  of  English  is  approxi- 
mately 1.3  bits  per  symbol,  thereby  indicating  an  achievable  data  com- 
pression on  English  text  of  3 to  1. 

d.  Two-Way  Channels  (T.  Cover,  R.  King,  A . El  Gamal , K . Jarett) 

An  attempt  is  being  made  to  determine  the  capacity  region 
of  the  Blackwell  two-way  channel. 
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B.  Project  6240.  OPTIMAL  CODES  WITH  VARIABLE  WEIGHT  SYMBOLS 

Principal  Investigator:  J.  Gill 

Staff:  N.  Cot 

1 . Objective 

The  objectives  of  this  project  are  to  determine  procedures  for 
producing  optimal  uniquely  decipherable  codes,  using  symbols  of  unequal 
costs,  and  to  study  the  properties  of  optimal  codes. 

2 . Current  Status  of  Work 

We  have  demonstrated  that  variable-length  coding  schemes  have 
applications  in  several  fields  other  than  communications.  These  include 

(a)  game  theory  where  they  are  associated  with  the  determina- 
tion of  optimal  gambling  strategies 

(b)  data  allocation  in  storage  and  data  retrieval  Cl 3 where 
their  application  is  based  on  the  relationship  between 
double-chained  trees  and  variable-length  codes  ; a file 
can  be  organized  as  a double -chained  tree  in  which  keys 
correspond  to  variable-length  code  words 

In  (b) , an  optimal  tree  (or  code)  corresponds  to  a minimal  av- 
erage search  time.  If  0 < b^  < b^  < . . . < b^  are  the  costs  of  the 
symbols  of  a code  alphabet  and  if  A is  determined  by 

*-bi  ,_bt  , 

A + A + • • • + A — 1 

then  the  converse  to  the  source -coding  theorem  provides  a lower  bound  on 
the  average  cost  c of  a code  [w^,...,wn)  with  code-word  probabilities 
(p^  > • • • «Pn  } : 

_ H(p^  , . • . , p^ ) 

C - log2  A 

where  H(p^,...,pn>  is  the  entropy  in  bits  of  the  probability  distribu- 
tion p.  , . . * ,p  * 

1 n 
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We  have  established  an  upper  bound  on  the  cost  of  optimal  codes 


H(p1,...,pn) 

C ' log2  A + b1  + 6 

where  & < b Is  a quantity  calculated  from  the  costs  bj,...,^.  Using 
this  upper  bound,  we  can  efficiently  find  nearly  optimal  codes  for  unequal 
probabilities . 

For  equiprobable  code  words,  a procedure  [2j  for  producing  opti- 
mal codes  has  been  analyzed  and  shown  to  require  OQog*'  n)  steps  to  pro- 
duce an  optimal  code  with  n code  words.  Current  research  deals  with 
extending  this  procedure  to  nearly  equiprobable  code  words. 

References 

1.  L.  Stanfel,  "Optimal  Trees  for  a Class  of  Information  Retrieval  Prob- 
lems," Info.  Stor.  Retr.,  9,  1973,  pp . 43-59. 

2.  N.  Cot,  "The  Complexity  of  the  Variable-Length  Encoding  Problem," 
Sixth  Southeastern  Conference  on  Combinatorics,  Graph  Theory,  and 
Computing,  1975. 

C.  Project  6302.  PROBLEMS  IN  CONTROL 

Principal  Investigator:  G.  F.  Franklin 

Staff:  G.  Aral 

1 . Objective 

This  study  is  concerned  with  design  problems  in  automatic  con- 
trol systems,  especially  when  model  uncertainty  is  present. 

2 . Current  Status  of  Work 

An  interpreter  in  APL  language  has  been  developed  for  linear 
control  design  problems.  Tills  system 

• allows  incorporation  of  new  routines  in  a trivial  way 
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• is  flexible  enough  to  accommodate  a wide  range  of  user 
familiarity  with  the  system,  allowing  use  of  its  func- 
tions in  either  a user-passive  or  user-active  manner 

This  transition  in  mode  is  not  abrupt  in  that  the  user  can  be  as  active 
as  his  knowledge  of  the  system  permits.  An  intensive  user  can  reduce 
terminal  time  as  he  becomes  more  familiar  with  the  system,  and  a one- 
time user  does  not  require  a detailed  study  of  function  descriptions. 

The  following  features  of  APL  make  it  attractive  as  an  inter- 
preter language . 

• It  is  designed  as  an  interactive  language. 

• This  language  has  powerful  algebraic  and  string-handling 
functions . 

• The  file  system  and  storage  capacity  are  extensive. 

The  three  types  of  functions  that  form  the  interpreter  are  design,  sys- 
tem, and  internal  functions.  Internal  functions  are  hidden  from  the 
user  and  are  virtually  inaccessible.  Design  and  system  functions  are 
niladic  and  act  on  predetermined  global  system  variables.  All  actions 
of  the  system  follow  three  steps: 

• possible  modification  of  the  content  of  the  globa’  system 
variables  to  be  used  by  the  functions 

• execution  of  specified  function(s) 

» possible  transfer  of  the  content  of  the  global  system 
variables  that  contain  the  results  of  the  executed  func- 
tion^) into  prescribed  global  user  variables 

All  the  commands  entered  by  the  user  are  analyzed  and  separated 
into  a function  name,  a separator,  and  a parameter  set.  The  separator  is 
either  a blank  or  dot;  for  example, 

SAMPL 

SAMPL  F,G, DELTA, PHI .GAMMA 
SAMPL  . 
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where  SAMPL  is  a function  name,  and  F, G, DELTA , PHI , GAMMA  is  the  parameter 
set . 

If  the  separator  is  a blank  and  the  parameter  set  is  empty,  the 
system  operates  in  the  passive  mode  and  guides  the  user  to  provide  the 
necessary  information  by  requesting  data  input  or  global  user  variable 
names  containing  the  data  and  requesting  the  output  medium  as  a global 
user  variable  name  or  display.  If  the  separator  is  blank  and  the  param- 
eter set  is  sufficient,  the  contents  of  the  given  input  variables  are 
transferred  to  the  global  system  variables,  the  function  is  executed, 
and  the  results  are  transferred  from  global  system  variables  to  the  cor- 
responding global  user  variables.  If  the  separator  is  a dot,  the  system 
executes  the  named  function  with  the  current  content  of  the  global  sys- 
tem variables  and  outputs  the  results  to  the  global  system  variables. 

The  following  features  of  the  system  are  a great  convenience. 

• The  user  can  set  up  and  execute  sequentially  a given  list 
of  commands. 

• A multilevel  activation  of  the  system  is  allowed  so  that, 
while  the  system  is  active,  a new  system  can  be  activated 
which  suspends  all  the  global  system  variables  of  the  ini- 
tial system  and  creates  a second  generation  of  global  sys- 
tem variables  . 

The  number  of  activated  levels  is  restricted  only  by  the  memory  of  the 
computer.  Return  to  the  upper  level  system  destroys  the  global  system 
variables  of  the  lower  level  system. 

Design  and  system  functions  are  incorporated  into  the  system. 
The  addition  of  many  other  functions  is  the  result  of  user  request  and 
programmer  effort. 


D.  Pro.ject  6418.  OPTICAL  COMPUTATION  BASED  ON  RESIDUE  ARITHMETIC 

Principal  Investigator:  J.  W.  Goodman 

Staff:  A.  Huang 


1 . Objective 

Since  1 October  1976,  this  project  has  been  reactivated.  Its 
new  goal  is  to  explore  the  feasibility  of  realizing  an  optical  computer 
based  on  residue  arithmetic. 
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Current  Status  of  Work 


Arithmetic  operations  based  on  a residue  number  system  have  the 
fundamental  advantage  of  not  requiring  carries  between  adjacent  columns. 
Operations  can  be  performed  in  parallel  channels,  with  no  need  for  one 
channel  to  wait  for  the  result  of  operations  performed  in  another  channel . 
The  inherent  parallelism  of  these  operations,  plus  the  well-known  capa- 
bility of  optical  systems  for  operating  in  many  parallel  channels,  has 
stimulated  interest  in  the  possibility  of  an  optical  computer  based  on 
residue  arithmetic  [11. 

If  the  value  of  a number  is  represented  by  the  spatial  location 
of  a small  spot  of  light,  the  operations  of  residue  conversion,  addition, 
subtraction,  multiplication,  and  return  to  a mixed  radix  system  can  be 
performed  by  certain  elements  called  "maps"  that  change  the  spatial  loca- 
tion of  the  spot  in  accord  with  the  operation  to  be  performed.  The  basic 
components  of  an  optical  residue  computer,  therefore,  are  prearranged  el- 
ements that  map  a set  of  possible  input  spot  locations  into  a set  of  out- 
put spot  locations  and  a device  that  selects  between  the  various  mapping 
elements.  To  add  two  numbers,  the  first  number  is  represented  by  the 
spatial  position  of  a light  spot,  and  a map  corresponding  to  addition  by 
the  second  number  is  selected.  The  position  of  the  output  light  spot 
represents  the  result  of  the  computation. 

An  arithmetic  unit  based  on  these  principles  has  certain  at- 
tractive properties  above  and  beyond  the  potential  speed  resulting  from 
the  parallelism  of  the  computations.  First,  the  accuracy  of  computations 
can  be  increased  simply  by  adding  new  modules  in  parallel;  second,  the 
computer  is  highly  fault  tolerant,  with  accuracy  degrading  slowly  with 
internal  defects  rather  than  failing  catastrophically. 

During  the  past  few  months,  we  have  examined  a variety  of  tech- 
nological approaches  to  realizing  the  maps  and  selectors;  however,  the 
ideas  for  organizing  and  structuring  the  mapping  operations  are  evolving 
too  raj:_.ily  for  us  to  select  one  particular  technological  approach.  Ac- 
cordingly, current  work  is  concentrated  on  improving  the  architecture 
envisioned  for  the  optical  computer. 

Spurting  in  February  1977,  this  work  will  be  supported  by  a 
contract  from  the  Air  Force  Office  of  Scientific  Research,  and  the  cur- 
rent JSEP  project  will  be  terminated. 
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Re ference 


1.  Alan  Huang,  "The  Implementation  of  a Residue  Arithmetic  Unit  via 

Optical  and  Other  Phenomena,"  Proc ■ of  the  1975  International  Opti- 
cal Computing  Conference,  Apr  1975. 


E.  Project  6503.  DATA  COMPRESSION  TECHNIQUES  AND  APPLICATIONS 

Principal  Investigator:  R.  M.  Gray 

Staff:  Y.  Linde 

1 . Objective 

Our  goals  are  to  develop  the  theory  and  computer-aided  design 
techniques  for  nearly  optimal  data-compression  algorithms. 

2 . Current  Status  of  Work 

With  the  joint  support  of  JSEP  and  APOSR,  a simulation  capa- 
bility for  trellis  coding  was  developed  on  the  Information  Systems  Lab- 
oratory Interdata  7/16  computer.  Numerous  simulations  were  conducted 
for  compressing  binary  and  gauss ian  sources  with  memory.  Decoders  for 
binary  sources  were  searched  exhaustively  up  to  length  four.  Decoders 
for  gaussian  sources  were  selected  via  ad  hoc  techniques,  particularly 
by  similarity  to  predictive  quantizer  decoders  and  from  Central  Limit 
Theorem  arguments.  The  gaussian  codes  performed  approximately  1 dB  bet- 
ter than  traditional  DPCM  techniques  and  suggested  some  additional  meth- 
ods that  should  produce  even  better  performance.  The  binary  codes  did 
not  yield  insight  into  good  code  design,  and  further  work  is  necessary. 
Because  the  Interdata  computer  is  now  gone,  some  time  will  be  required 
to  adopt  the  programs  to  the  new  computer,  and  then  the  simulation  stud- 
ies can  be  continued.  In  the  meantime,  our  focus  will  be  directed  to- 
ward analysis  of  the  newly  developed  coding  structures. 
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F.  Project  6601.  INFORMATION  THEORY  AND  INFORMATION  PROCESSING 

Principal  Investigator:  M.  E.  Heilman 

Staff:  L.  Pinto  de  Carvalho,  R.  Kahn 

1 . Objective 

The  objective  of  this  project  is  to  investigate  the  relation- 
ships between  information  theory,  information  processing,  and  complexity 
theory . 


2.  Current  Status  of  Work 

Coding  schemes  for  the  wiretap  channel  with  feedback  are  being 
considered.  A trade-off  between  entropy  production  (to  confuse  the  wire- 
tapper)  and  information  flow  has  been  established.  This  trade-off  and 
recent  results  from  other  multiuser  communication  channel  models  are 
being  studied  for  use  in  bounding  the  achievable  rate-equivocation  region. 


G.  Project  6701.  ADAPTIVE  FILTERING 

Principal  Investigator:  B.  Widrow 

Staff:  C.  Williams,  M.  Larimore,  J.  Treichler 

1 . Objective 

The  objective  of  this  project  is  to  study  the  effectiveness  of 
various  adaptive  techniques  in  digital  communication  systems.  Interest 
has  focused  on  adaptively  computed  direct-channel  models  in  systems  where 
doppler  effects,  multipath,  and  additive  noise  are  present.  Other  work 
is  concerned  with  the  application  of  adaptive  techniques  to  the  detection 
of  spectral  lines  in  noise. 

2 . Current  Status  of  Work 

The  study  of  adaptive  channel  modeling  in  the  presence  of  dop- 
pler and/or  clock  asynchronizat ion  is  continuing  in  a new  direction.  We 
have  obtained  a good  characterization  of  the  LMS-driven  adaptive  model 
performance  and  are  now  investigating  schemes  for  augmenting  the  algorithm 
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to  Improve  its  asynchronized  performance.  These  schemes  are  divided  into 
(a)  adaptive  techniques  for  adjusting  the  receiver  clock  based  on  model 
behavior  and  (b)  techniques  for  augmenting  the  algorithm  to  include  an 
adaptive  correction  in  the  channel  model.  This  line  of  research  is  new 
and  appears  to  be  promising. 

The  extraction  of  a narrowband  signal  buried  in  noise  is  one 
of  the  oldest,  yet  still  important,  problems  in  the  areas  of  communica- 
tion and  signal  processing.  We  are  continuing  our  development  of  the 
adaptive  line  enhancer  (ALE)  to  separate  a signal  into  its  narrow-  and 
broadband  components.  A conclusive  study  of  the  line  enhancer,  using  an 
all-zero  (feedforward)  adaptive  filter,  has  been  completed.  We  are  now 
investigating  a line  enhancer  that  employs  a feedback  and  feedforward 
("poles"  and  "zeros")  adaptive  filter.  The  poles  provide  a high  "q" 
filter  with  a small  number  of  adaptive  weights  and,  thereby,  is  highly 
suited  for  the  detection  of  narrowband  signals  buried  in  noise.  The  op- 
eration of  the  adaptive  feedback  filter  is  significantly  more  involved 
than  that  of  the  feedforward  filter  and  is,  of  itself,  a research  topic. 
Our  preliminary  results  indicate  that  the  feedback  line  enhancer  is  a 
great  improvement  over  the  feedforward  device. 


H.  Project  6851.  IMAGING  OF  GAMMA-RAY  DISTRIBUTIONS 

Principal  Investigator:  A.  Macovski 

Staff:  A.  L.  Steinbach,  D.  Rosenfeld 

1 . Objective 

The  objectives  of  this  work  are  to  study  the  fundamental  char- 
acteristics of  a gamma -ray  source  distribution  and  to  derive  optimal 
methods  for  imaging  such  distributions.  Because  refractive  effects  do 
not  exist  at  these  energies,  various  attenuative  masking  structures  are 
being  investigated  with  respect  to  efficiency , resolution , depth  deline- 
ation, and  s ignal-to-noise  ratio. 

2 . Current  Status  of  Work 

This  represents  the  conclusion  of  work  on  this  project.  Two 
doctoral  dissertations  are  being  prepared.  Tbey  will  describe  the  use 
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of  zone  plates  and  modulated  apertures  to  efficiently  image  gamma-ray 
sources . 

In  the  work  based  on  zone  plates,  the  three-dimensional  point 
response  of  the  entire  system  has  been  analyzed,  including  encoding  and 
decoding.  The  output  of  the  computer  plotter  graphically  demonstrates 
system  performance,  including  lateral  and  depth  resolution  and  sidelobe 
response.  Several  zone-plate  configurations  have  been  studied  to  deter- 
mine those  that  provide  the  best  performance  in  various  dimensions.  Ap- 
odization  has  been  investigated  as  a trade-off  between  the  size  of  the 
central  lobe  and  the  amplitude  of  the  sidelobes.  The  results  provide  a 
wealth  of  design  criteria  for  these  structures. 

In  the  work  using  modulated  apertures,  many  computer  simula- 
tions and  analytical  studies  have  been  performed.  The  analyses  have 
demonstrated  that  the  resultant  reconstruction  is  optimum  in  the  maximum 
likelihood  sense.  A method  has  been  derived  to  determine  the  optimal 
code  parameters  for  a given  ob/ject  size;  this  code  parameter  represents 
the  number  of  imaging  apertures  open  at  any  given  time  period  of  the 
code  cycle.  Several  codes  have  been  considered  and  implemented.  A num- 
ber of  criteria  have  been  used,  including  maximizing  the  maximum  signal- 
to-noise  ratio  in  the  image  or  maximizing  the  minimum  signal-to-noise 
ratio.  A simulated  reconstruction  of  a three-dimensional  object,  in  the 
presence  of  counting  noise,  is  under  way  to  indicate  the  ability  to  sep- 
arate planes  in  space. 


I.  Project  7050.  STUDIES  IN  STATISTICAL  SYSTEM  THEORY 

Principal  Investigators:  T.  Kailath,  M.  Morf 

Staff:  J.  Dobbins,  B.  bevy,  S-Y.  Kung , E.  Verriest , 

A . Vie ira 


1 . Objective 


The  purpose  of  this  project  is  to  study  basic  problems  in  sta- 
tistical system  theory. 
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Current  Status  of  Work 


2 . 

a .  Square-Root  Algorithms  for  the  Continuous-Time  Linear 
Least-Squares  Estimation  Problem  (M.  Morf,  B.  Levy, 

T.  Kailath) 

We  have  found  a simple  differential  equation  to  propagate 
the  triangular  square  root  of  the  error  covariance  of  the  linear  estimate 
of  the  state.  This  differential  equation  does  not  involve  the  antisym- 
metric matrix  which  appeared  in  the  previous  work  of  Andrews  and  Tapley 
and  Choe . We  combine  both  results  by  constructing  an  explicit  antisym- 
metric matrix  that  generates  the  orthogonal  transform  that  continuously 
triangularizes  the  square-root  factor. 

As  in  the  discrete  case,  an  array  method  for  handling  the 
continuous  square  root  was  developed  and  yields  the  usual  norma  1 Nations 
of  stochastic  calculus  and  the  results  of  the  other  approaches.  In  the 
constant  model  case,  it  is  demonstrated  that  the  Chandrasekhar  equations 
are  a set  of  square-root  equations  for  the  derivative  of  the  error  vari- 
ance . 


b.  Multichannel  Maximum  Entropy  (M.  Morf,  A.  Vieira, 

T.  Kailath) 

Our  work  on  multichannel  maximum  entropy  methods  (MMEM)  is 
progressing.  A paper  has  been  completed  and  has  attracted  considerable 
interest.  Pusey  (at  EXON)  has  implemented  our  algorithm  in  the  computer 
and  verified  the  expected  well-behavior  of  the  numerical  properties.  The 
derived  spectral  estimates  appear  to  be  very  close  to  the  maximum  like- 
lihood estimates  (quote  from  F.  Muir,  EXON,  a frequent  guest  in  the  Stan- 
ford Geophysics  Department).  An  invitation  has  been  extended  to  Profes- 
sor M.  Morf  to  present  a paper  on  this  subject  at  the  1977  JACC  Conference 
in  San  Francisco  to  a special  topics  session  on  geophysical  data  proces- 
sing organized  by  J.  Mendel  (USC) . 

c.  Ladder  Structure  for  Least-Squares  Estimation 
(M.  Morf,  A.  Vieira,  T.  Kailath) 

Ladder  structures  increasingly  appear  to  be  the  "natural" 
canonical  form  for  estimation  problems.  We  have  developed  exact  least- 


11 


squares  recursions  for  such  relations  (including  multidimensional  ver- 
sions) and  have  found  that  they  have  even  less  computations  and  storage 
requirements  than  our  earlier  forms  of  fast  algorithms  for  linear  pre- 
diction . 

These  recursions  do  not  require  "auxiliary"  variables  as 
did  the  earlier  forms,  except  for  a set  of  scalars  that  perform  the 
function  of  "optimal  weights"  in  real-time  correlators  that  compute  the 
model  parameters.  In  a project  on  speech  processing,  we  have  written  a 
survey  paper  in  which  some  of  these  new  results  are  included.  They  will 
be  extremely  useful  in  that  project  because  they  will  be  part  of  a set 
of  fast  and  efficient  speech-encod ing  algorithms  for  high  quality  and 
low  data-rate  properties.  A more  extended  version  of  this  paper  on 
scalar  ladder  forms  is  planned  for  the  IEEE-ASSP  Journal  and  a multi- 
channel version  for  the  IEEE-AC  Journal. 

d . Square-Root  Algorithms  for  Model  Sensitivity  Analysis 

(M.  Morf,  T.  Kailath,  J.  Dobbins,  E.  Verriest) 

We  have  developed  square-root  algorithms  for  calculating 
the  true  error  covariance  of  a Kalman  filter  estimate  based  on  an  in- 
correct model  of  the  process  under  observation.  These  algorithms  have 
been  extended  to  the  smoothing  problem.  When  based  on  an  incorrect 
model,  the  filter  arrays  do  not  properly  orthonorma 1 ize  the  estimate 
residuals.  Our  algorithms  propagate  on  additional  arrays  that  complete 
the  orthonormalization. 

We  are  applying  these  methods  to  the  singular  perturba- 
tion problem.  Work  is  also  progressing  on  connecting  the  successful 
square-root  and  scattering-theory  frameworks  for  least-squares  estima- 
tion . 

e . Minimal  Partial  Realizations  of  M I MO  Systems 

(S-Y.  Kung,  T.  Kailath,  M.  Morf) 

We  have  demonstrated  that  a modified  Lanczos  recursion 
can  be  used  to  attain  a minimal  partial  realization  of  SISO  linear  sys- 
tems. The  first  step  is  to  define  a Hankel  matrix  H = th  ^1  where 

th  i+J 

h^  is  the  (i+j)  Markov  parameter.  The  key  step  is  to  transfer  this 
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Hankel  matrix  into  an  echelon  matrix  until  a dependent  column  is  obtained. 
The  denominator  of  the  solution  is  then  solved  by  noting  that  its  coeffi- 
cients are  simply  those  dependent  coefficients  of  the  dependent  column. 
This  procedure  guarantees  that  the  order  of  our  solution  system  (say  n) 
is  indeed  a minimal  one. 

We  will  extend  the  above  recursive  algorithm  to  the  MIMO 
partial  realization  problems.  A pseudo-canonical  polynomial  matrix  is 
first  introduced.  This  Popov  form  has  an  indispensable  structural  prop- 
erty and  very  powerful  applications  in  many  design  problems  in  addition 
to  the  partia 1-rea 1 izat ion  one.  Our  procedure  is  again  to  transform  the 
Hankel  matrix  into  an  echelon  form  and  to  find  m (where  m is  the  num- 
ber of  inputs)  primary  dependent  columns.  It  can  be  shown  that  this  pro- 
cedure guarantees  a set  of  minimal  input  dynamical  indices  (controllabil- 
ity indices)  and,  accordingly,  a minimal  order  of  the  solution  system. 

f . New  Results  on  Two-Dimensional  Systems  Theory 

(B.  Levy,  S-Y.  Kung , M.  Morf,  T.  Kailath) 

In  recent  years,  linear  system  theory  has  been  applied 
intensively  to  estimation  and  control  problems.  At  the  same  time,  image 
processing  has  attracted  increasing  interest.  As  a result,  attempts  have 
been  made  to  extend  the  techniques  of  systems  theory  to  two-dimensional 
problems  either  from  a polynomial  ( input -output ) or  state-space  (inter- 
nal) point  of  view. 

In  the  first  part  of  our  work,  we  examined  some  new  re- 
sults obtained  from  two-dimensional  polynomial  matrices  and  have  suc- 
cessfully extended  the  existing  one-dimensional  results  on  GCRD  extrac- 
tion, Sylvester  resultants,  and  matrix  fraction  descriptions  (MFD)  to 
the  2-D  matrices.  Some  of  our  results  appear  to  be  unique  for  such  mul- 
tidimensional problems  as  the  existence  of  "primitive  factorizations"  in 
addition  to  general  factorizations  of  the  2-D  polynomial  matrix. 

In  the  second  part  of  our  work,  we  examined  some  defini- 
tions of  state  controllability  and  observability  and  their  relationship 
to  the  minimality  of  2-D  systems.  We  have  also  investigated  new  circuit- 
realization  methods,  using  a minimal  number  of  dynamic  elements  and  a 
new  technique  for  2-D  digital-filter  hardware  implementations  of  2-D 
transfer  functions. 
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We  have  also  examined  the  2-U  systems  from  a computer 
algorithm  point  of  view.  Our  results  include  a generalization  of  Lev- 
inson's algorithm,  finding  a Hankel  matrix,  and  developing  a recursive 
realization  algorithm  given  a set  of  Markov  parameters.  There  are  still 
many  open  problems  for  2-D  systems,  such  as  stability  tests,  innovations 
representation,  and  Kalman  filtering. 
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II.  DIGITAL  SYSTEMS 


A.  Project  6961.  DESCRIPTION  LANGUAGES  AND  DESIGN  FOR  GENERAL  PURPOSE 

COMPUTER  ARCHITECTURES 

Principal  Investigators:  M.  J.  Flynn, 

W.  M.  vanCleemput 

Staff:  T.  Bennett,  D.  Hanson,  J.  Hupp,  R.  Lee, 

K.  Stevens,  S.  Wakefield 

1 . Objective 

The  purposes  of  this  study  are  to  develop  description  languages 
to  describe  computer  architectures  and  to  establish  a basis  for  under- 
standing the  limits  of  computer-architecture  design  (to  determine  the 
fastest  speed  for  the  execution  of  a program)  . 

2 . Current  Status  of  Work 

a . Evaluation  of  Existing  Hardware  Description  Languages 
(D.  Hanson,  W.  vanCleemput) 

To  evaluate  the  feasibility  of  describing  computer  archi- 
tectures by  means  of  a hardware  description  language,  several  of  these 
languages  were  studied.  From  a large  number,  three  were  selected  based 
on  the  availability  of  a working  compiler  and  simulator  and  on  the  suit- 
ability of  the  language  for  multilevel  description  of  digital  systems. 
Compilers  for  the  CDL  (Chu,  University  of  Maryland)  and  the  CASSANDRE 
(Mermet , University  of  Grenoble,  France)  were  implemented,  and  these 
languages  were  studied  and  compared  in  terms  of  their  descriptive  power. 

A compiler  for  the  DDL  (Dietmeyer,  University  of  Wisconsin)  was  obtained 
and  is  being  adapted  for  the  IBM  370  system. 

b.  Development  of  a Structural  Description  Language 
TwT  vanCleemput) 

Most  of  the  existing  hardware  description  languages  em- 
phasizes a behavioral  description  and  neglects  the  structural  aspects  of 
a design.  A behavioral  description  is  often  more  than  sufficient  to  de- 
scribe an  existing  architecture;  however,  to  aid  in  the  design  process, 
it  is  necessary  to  capture  all  available  information — structure  as  well 
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as  (intended)  behavior.  The  SDL  language  was  developed  to  describe  ac- 
curately the  structural  properties  of  the  system,  and  a compiler  is  cur- 
rently operational.  A special  characteristic  of  the  SDL  is  that  it  is 
hierarch ial  in  nature. 

c . Applications  of  the  Structural  Design  Language 
(T.  Bennett,  J.  Hupp,  K.  Stevens,  W.  vanCleemput) 

One  application  of  a structural  design  language  is  to 
serve  as  the  input  medium  for  a physical  design  system.  Such  a system 
may  consist  of  several  subsystems  to  perform  such  tasks  as  printed-cir- 
cuit layout,  integrated-circuit  layout,  gate-level  logic  simulation, 
circuit-level  simulation,  fault  test  generation,  and  automated  logic 
diagram  generation.  A printed-circuit  layout  system  is  being  imple- 
mented, and  a prototype  will  soon  be  operational . An  integrated-circuit 
layout  system  is  in  the  design  phase.  Interfaces  between  the  SDL  lan- 
guage and  several  existing  logic-  and  circuit-level  simulators  are  in 
the  planning  phase. 

d.  Bounds  for  Maximal  Parallelism  (R.  Lee,  M.  Flynn) 

This  is  a study  of  the  performance  limits  of  a single 
program  executed  on  a large  number  of  identical  processors  operating  in 
parallel  in  an  MIMD  (multiple  instruction  multiple  data)  organization. 
With  the  rapidly  decreasing  cost  of  LSI  microprocessors,  it  is  now  eco- 
nomically feasible  to  consider  an  army  of  processors  within  the  computer 
architecture  to  speedup  a computation,  even  at  reduced  efficiency  of 
each  component  processor.  No  longer  i<  the  processor  the  hallowed  CPU 
(central  processing  unit)  or  the  most  valuable  resource  to  be  utilized 
with  the  greatest  efficiency.  Some  "acceptable"  level  of  efficiency 
should  be  obtained,  however,  but  we  must  determine  what  type  of  speedup 
can  be  expected  by  increasing  the  number  of  processors  even  if  we  ignore 
the  problems  of  control  and  communication  that  accompany  the  cooperation 
and  competition  between  these  processors. 

We  first  defined  a general  model  computation  on  a p-par- 
allel  processor,  distinguishing  clearly  between  the  logical  parallelism 
(p*  processors)  inherent  in  a computation  and  the  physical  parallelism 
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(p  processors)  available  in  the  computer  organization.  We  then  deter- 
mined such  performance  measures  as  execution  time,  speedup,  efficiency, 
and  space-time  product  by  which  we  could  evaluate  the  performance  im- 
provements (if  any)  of  p-parallel  processor  systems  over  uniprocessor 
systems.  The  results  obtained  indicated  that,  generally,  performance 
depends  on  both  the  computer  architecture  and  the  computation. 

We  then  derived  necessary  and  sufficient  conditions  for 
the  maximum  attainable  speedup  of  a p-parallel  processor  over  a unipro- 
cessor to  be  S < min  (p/ln  p,  p*/ln  p*).  Despite  the  many  different 
P 

vie w s concerning  the  potential  speedup  of  parallel  processor  systems, 
this  bound  has  never  before  been  established. 

In  addition,  if  there  are  always  sufficient  processors 
(p  > p*),  conditions  can  be  derived  under  which  the  maximum  attainable 
speedup  of  a computation  is  p*/ln  p* , maximum  efficiency  is  1/ln  p* , 
minimum  execution  time  is  X In  p*/p*,  and  minimum  space-time  pro- 
duct is  X In  p*  (where  T is  the  execution  time  of  the  computa- 

tion on  a uniprocessor).  In  fact,  the  empirical  speedups  obtained  for 
a large  number  of  different  computations  indicated  that  80  percent  of 
all  programs  examined  satisfied  these  conditions  and  had  maximum  speed- 
ups  of  less  than  p*/ln  p* . 

e . Parallel  Information  Processing  in  Biological  Systems 
(S.  Wakefield,  M.  Flynn) 

The  interconnections  and  types  of  synapses  between  units 
of  a particular  neural  subsystem  (the  stomatogastric  ganglion  of  the 
lobster)  have  been  determined  by  biologists,  as  have  the  stereotyped 
motor  patterns  that  it  produces;  however,  the  exact  mechanisms  and  se- 
quence and  duration  bounds  of  impulse  bursts  that  must  under 1 ie  the  pro- 
duction of  the  coordinated  muscle-activation  patterns  are  unknown.  Such 
a mechanism  would  be  analogous  to  the  switching  network  responsible  for 
the  traversal  of  states  in  a digital  sequential  circuit.  Because  of 
this  analogy,  this  and  similar  simple  biological  information-processing 
subsystems  are  being  investigated.  In  addition,  the  speed,  power  re- 
quirements, size,  information  capacity,  and  other  characteristics  of 
individual  neurons  are  being  compared  to  those  of  electronic  informa- 
tion-processing components. 
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B.  Project  7151.  COMPUTER  ARCHITECTURE 

Principal  Investigators:  E.  J.  McCluskey, 

S.  S.  Owicki 

Staff:  J.  Losq,  R.  C.  Ogus  , D.  J.  Rossetti, 

A.  Spector,  F.  W.  Terman 

1 . Objective 

The  objective  of  this  study  is  to  gain  new  insight  into  com- 
puter architecture  and  operating  systems.  We  intend  to  design  a small 
operating  system  and  to  verify  its  correctness  . Another  goal  is  to  de- 
velop new  methods  to  evaluate  system  performance  and  to  use  the  results 
obtained  for  better  system  design. 

2 . Current  Status  of  Work 

This  project  has  developed  in  two  major  directions — an  inves- 
tigation of  operating-system  correctness  and  the  development  of  analyt- 
ical models  and  monitoring  systems  for  performance  evaluation. 

The  bas,c  methods  for  designing  and  verifying  the  correctness 
of  concurrent  programs  have  been  analyzed.  The  design  of  a probably 
correct  operating  system  is  under  way,  and  the  high-level  description 
is  nearing  completion. 

A set  of  trace  programs  has  been  developed  to  measure  micro- 
processor performance  and  to  provide  unique  information  concerning  the 
programming  ease  of  various  microprocessors.  The  tracing  is  highly 
powerful  (can  record  interrupts  and  I/O  behavior)  and  almost  transparent 
(external  hardware  reduces  significantly  the  amount  of  perturbation  in- 
troduced by  the  tracer)  . 

The  investigation  of  interleaved  memory  is  being  extended  to 
the  effects  of  request  ordering  and  to  the  most  efficient  ordering.  A 
new  model  has  been  developed  to  evaluate  the  performance  of  gracefully 
degradable  multiprocessor  systems,  based  on  such  measures  as  availabil- 
ity, processing  power,  and  proportion  of  time  in  the  degraded  mode.  It 
also  allows  one  to  determine  the  best  trade-off  between  computing  p.  wer 
and  availability  or  between  hardware  vs  software  fault  detection. 

A new  testing  method  indicates  that  random  inputs  can  be  used 
efficiently  to  test  sequential  circuits.  A statistical  analysis  revealed 
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that  this  method  is  potentially  the  most  cost  effective  for  in-the-field 
testing . 


a.  Specification  and  Verification  of  Monitors  (S . Owicki) 

A monitor  is  a programming  language  construction  that  de- 
fines a logically  related  group  of  shared  data  items  and  a set  of  opera- 
tions on  those  items.  The  operations  are  the  only  means  by  which  programs 
may  access  the  shared  data,  and  the  monitor  includes  synchronization  to 
ensure  that  processes  do  not  interfere  with  each  other  as  they  perform 
their  operations.  This  limited  access  simplifies  the  verification  task. 

A monitor  is  specified  by  listing  the  shared  data  items, 
their  initial  values,  and  the  effects  of  each  operation.  Verification  of 
programs  is  simplified  by  describing  the  operations  in  terms  of  variables 
that  are  private  to  the  program  invoking  the  operation.  The  relationship 
between  the  private  and  shared  variables  is  expressed  by  an  invariant  re- 
lation, which  is  true  for  the  initial  monitor  values  and  is  preserved  by 
each  operation. 

There  are  two  steps  in  the  process  of  verifying  systems 
that  use  monitors.  The  first  shows  that  the  monitor  satisfies  its  spec- 
ifications (that  the  operations  preserve  the  invariant  and  have  the  re- 
quired effect  on  the  private  variables),  and  the  second  verifies  the  pro- 
cesses that  call  the  monitor  by  using  the  specifications  of  the  monitor 
onerations.  In  some  instances,  the  correct  behavior  of  one  process  de- 
pends on  the  actions  of  another.  Here  it  is  convenient  to  use  a process 
invariant — an  invariant  relation  on  the  private  process  variables  . These 
tools  appear  to  be  powerful  enough  for  most  applications  of  concurrent 
programming  with  monitors.  Further  work  is  required  to  develop  methods 
for  treating  dynamically  allocated  resources  that  do  not  fit  the  monitor 
pattern . 


b.  Operating  System  Design  (A.  Spector) 

Our  work  has  focused  on  the  design  of  a paging  interactive 
time-sharing  system  suitable  for  a computer  such  as  a PDP11/70.  Although 
we  do  not  expect  to  include  very  complex  I/O  or  protection  mechanisms,  we 
are  attempting  to  develop  a realistically  usable  system. 
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Tho  system  will  be  a hierarchy  consisting  of  two  levels — 
the  kernel  and  the  supei'visor.  The  kernel  is  the  more  primitive  run- 
time support  system  that  implements  the  basic  operations  of  the  high- 
level  programming  language  used  by  the  supervisor.  More  specifically, 
it  contains  the  machine-dependent  operating-system  independent  procedures 
that  provide  concurrent  processes  and  monitors.  It  also  translates  cer- 
tain privileged  hardware  operations  to  language-level  primitives.  For 
example,  because  supervisor  processes  use  monitors  to  communicate,  the 
kernel  must  maintain  the  queues  necessary  to  support  critical  sections 
and  wait  and  signal  primitives.  An  unusual  feature  is  that  the  kernel 
maintains  detailed  per-process  state  information  that  can  be  used  by 
such  supervisor  routines  as  scheduling  and  accounting. 

Although  the  nucleus  is  small  and  can  be  thought  of  as  a 
single  entity,  the  supervisor  performs  many  complex  functions  and  is  best 
considered  as  a hierarchy.  For  example,  the  long-term  process  scheduler 
is  a primitive  operation  of  the  supervisor  and  underlies  the  operation 
of  the  whole  system.  Our  design  calls  for  this  process  to  run  synchro- 
nously and  to  communicate  periodically  to  the  kernel  a set  of  processes 
that  can  be  executed  in  some  ensuing  interval.  The  memory  allocation 
process  (also  one  of  the  basic  operations)  contains  the  logic  to  allocate 
physical  memory  to  processes  and  to  support  shared  and  r.onshared  pages. 
Scheduling  and  memory-allocation  decisions  can  be  made  with  the  use  of 
state  information  recorded  by  the  kernel  and  from  information  contained 
in  certain  supervisor  monitors.  Higher  levels  of  the  supervisor  hier- 
archy include  sections  to  support  I/O  and  user  program  calls. 

During  this  reporting  period,  our  work  has  focused  on  the 
design  of  the  most  primitive  portions  of  the  system.  We  believe  that 
such  issues  as  the  accurate  assignment  of  meanings  to  language  primitives 
and  the  proper  definition  of  the  boundary  between  the  kernel  and  super- 
visor will  greatly  influence  the  outcome  of  our  project.  As  a result, 
we  are  proceeding  with  caution  in  the  hope  of  preparing  a sound  basis  on 
which  to  construct  a carefully  structured  and  verified  system. 

c.  Trace  Facility  (D.  J.  Rossetti) 

The  STRAP/370  instruction  tracing  facility  described  in 
previous  reports  is  providing  data  for  various  studies  in  the  areas  of 
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performance,  architecture,  and  operating  systems.  A paper  entitled  "The 
Design  and  Implementation  of  an  Operating  System  Tracer"  has  been  sub- 
mitted for  publication  in  the  Communications  of  the  Association  for  Com- 
puting  Machinery.  It  describes  the  unique  aspects  of  the  tracer,  pro- 
vides some  examples  of  its  use,  and  recommends  further  applications  of 
the  data  and  the  technique. 

Our  approach  is  being  extended  into  the  microprocessor 
area,  retaining  most  of  the  desirable  features  of  the  original  technique. 
The  purpose  is  to  evaluate  computer  architectures  from  the  standpoint  of 
development  and  programming  ease  as  opposed  to  processing  speed  because 
software  development  has  become  a major  cost  in  microprocessor  use.  The 
principal  tool  is  a detailed  level  instruction  tracer  having  the  follow- 
ing characteristics. 

• It  is  independent  of  the  program  being  measured  to  the 
extent  that  one  microprocessor  can  be  used  to  measure 
another . 

• The  instruction  stream  is  sampled  in  bursts  at  a rate 
determined  by  the  user. 

• Minimal  perturbation  is  introduced  into  the  measured 
system  because  the  trace  is  collected  with  hardware 
assistance  and  not  by  interpretation. 

• Actual  systems  and  programs  can  be  measured,  including 
I/O  and  interrupts ; the  trace  is  not  gathered  in  an 
artificial  environment,  such  as  a simulator. 

The  tracer  has  been  used  to  gather  instruction  frequency 
data  for  application  programs  in  our  microprocessor  laboratory.  Current 
plans  are  to  collect  extensive  trace  data  for  studies  in  architecture 
evaluation  and  to  gain  understanding  of  microprocessor  program  and  in- 
put/output behavior. 

d . Memory  Interleaving  ( F . W . Terma n ) 

A series  of  models  of  interleaved  memory  systems  has  been 
investigated  by  means  of  a trace-driven  simulation.  The  basic  model  ex- 
tends the  one  developed  by  Burnett  and  Coffman  [1 1 to  the  architecture 
of  the  IBM  360/370  with  its  variable-length  instructions  and  operands  ; 
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it  also  includes  multibyte  transfers  per  memory  access.  Variations  of 
the  basic  simulation  model  facilitates  study  of  the  effects  of  channel 
interference,  data  request  reordering,  and  data  queue  emptying. 

Memory  requests  for  the  simulation  are  obtained  from  two 
sets  of  instruction-by-instruction  trace  records.  The  first  set,  pro- 
duced by  the  program  TRACE/360  written  by  Terman  , traces  the  problem 
state  component  of  typical  programs  running  on  the  IBM  360/370.  The 
second  set,  produced  by  the  program  STRAP/370  written  by  Rossetti, 
traces  samples  of  the  total  activity  of  the  CPU,  including  both  super- 
visor and  problem  states  . 

The  theoretical  predictions  of  Burnett  and  Coffman  for 
the  increase  in  memory  bandwidth  caused  by  interleaving  are  found  to 
fit  well  with  the  simulation  results  obtained  for  the  fetching  of  in- 
structions. The  usefulness  of  fetching  instructions  in  blocks  is  lim- 
ited by  the  relatively  high  frequency  of  branches  on  the  IBM  360/370 
machines.  Consequently,  an  active  channel  has  relatively  little  effect 
(less  than  10  percent  degradation)  on  the  fetching  of  instructions  from 
an  interleaved  memory . 

For  the  transfer  of  operands  to  and  from  memory,  the 
simulation  results  reveal  only  one-half  the  increase  in  memory  bandwidth 
predicted  by  the  analysis  of  Burnett  and  Coffman.  This  indicates  that 
data  references  on  the  IBM  360/370  are  not  random  as  were  assumed  . The 
effect  of  an  active  channel  is  again  small  because  of  the  number  of  idle 
modules  in  a typical  memory  cycle.  A larger  degradation  (10  to  20  per- 
cent) occurs  if  the  data  requests  are  forced  to  "catch  up"  with  the  in- 
struction requests  before  a successful  branch  is  executed.  This  is  a 
reasonable  restriction  because  the  branch  address  cannot  be  calculated 
with  certainty  until  the  values  of  the  general-purpose  registers  are 
known.  On  the  other  hand,  a significant  improvement  (15  to  30  percent) 
in  the  memory  bandwidth  can  be  obtained  by  allowing  the  data  requests  to 
be  filled  out  of  order.  These  effects  are  being  analyzed  in  detail  and 
compared  to  the  predictions  of  other  theoretical  models. 

This  investigation  is  being  extended  to  include  the  ef- 
fects of  interference  between  multiple  CPUs.  The  effects  of  the  inter- 
dependencies between  instructions  and  operands  and  between  the  operands 
themselves  will  also  be  examined. 
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e.  Gracefully  Degradable  Systems  and  Their  Performance 
(.J.  Losq) 

The  recent  development  of  multiprocessor  systems  that 
offer  tolerance  to  failures  by  switching  to  a degraded  mode  of  operation 
following  failure  occurrences  has  the  vast  potential  of  providing  large 
processing  power  coupled  with  high  reliability  (availability).  A general 
model  has  been  developed  T2l  to  evaluate  the  performance  of  such  ; stems. 
It  takes  into  account  the  architecture  of  the  hardware,  characteristics 
of  the  various  fault  detection  and  recovery  mechanisms,  unreliability  of 
the  hardware  components,  and  failure  rate  of  the  operating  system.  The 
model  provides  such  measures  of  overall  system  performance  as  average 
processing  power,  mean  time  between  crashes,  and  the  proportion  of  time 
spent  by  the  system  in  the  degraded  mode.  The  best  trade-offs  between 
high  processing  power  and  availability  (or  between  hardware  vs  software 
fault  detection)  are  obtained. 

This  model  can  also  determine  overall  performance  when 
computer  systems  are  seen  as  a combination  of  hardware  and  software,  and 
it  can  be  applied  to  almost  all  current  multiprocessors.  A detailed  ex- 
ample emphasized  the  large  contribution  of  software  errors  to  the  system 
down  time.  Future  work  will  concentrate  on  an  analysis  of  a recovery 
strategy  and  on  the  general  problem  of  fault-tolerant  software. 

f . Efficiency  of  Compact  Testing  for  Sequential  Circuits 

(J.  Losq) 

Compact  testing  uses  random  inputs  to  test  digital  cir- 
cuits. Failures  can  be  detected  by  comparison  between  some  statistical 
property  of  the  circuit  (for  example,  the  frequency  of  logic  l's  on  the 
output  lines)  and  the  same  property  for  the  fault-free  circuit.  Although 
such  a method  is  commonly  used  in  the  field,  no  previous  study  indicates 
its  usefulness  and  limitations. 

This  investigation  f3  ] revealed  that  compact  testing  is 
indeed  efficient  although  it  has  some  inherent  limitations.  Most  of  the 
failures  in  a sequential  circuit  can  be  detected;  however,  compact  test- 
ing cannot  fully  guarantee  that  a circuit  is  fault  free.  Synchronization 
is  not  a major  problem — any  long  sequence  of  random  inputs  will  act  as  a 
synchronizing  sequence,  and  most  circuits  can  be  synchronized  in  such  a 
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manner  in  a matter  of  seconds.  The  great  majority  of  failures  inside 
the  memory  elements  can  be  detected,  even  with  short  tests.  Compact 
testing  also  detects  most  of  the  failures  in  the  combinational  parts; 
its  efficiency  is  largely  dependent  on  the  test  length  and  on  the  char- 
acteristics of  the  random-number  generators. 

g.  Monitors  for  Signal  Activity  (R.  Ogus) 

Probabilistic  models  have  been  developed  [41  to  describe 
logic  circuits  where  a probability  can  be  assigned  to  a signal  line  to 
indicate  that  the  signal  is  a logic  "l"  given  the  probability  of  the  in- 
puts of  the  circuit  being  a ”l.”  Algorithms  have  been  devised  to  derive 
the  signal-probability  expressions  as  functions  of  the  input-signal  prob- 
abilities. The  signal  probability  can  be  used  in  computing  the  probabil- 
ity of  detecting  faulty  behavior  in  circuits  and  also  in  the  calculation 
of  the  signal  reliability  of  a circuit  (the  probability  that  the  circuit 
output  is  correct).  In  addition,  a scheme  has  been  described  to  select 
efficient  test  sets  from  random  inputs  [51 . This  method  measures  the 
effect  on  a circuit  output  of  exercising  any  particular  input  variable, 
and  the  inputs  can  be  weighted  according  to  their  importance.  This  re- 
quires measuring  both  the  input-  and  output-signal  probabilities.  Other 
uses  of  the  measure  are  found  in  error  latency  studies  and  in  system 
testing . 

In  many  cases,  circuits  are  too  large  to  be  studied  ana- 
lytically and,  as  a result,  it  is  desirable  to  have  a hardware  monitor 
unit  that  can  actually  measure  signal  activity  and  display  the  signal 
probability.  Two  monitor  designs  have  been  proposed.  The  first  is  a 
simple  stand-alone  single-probe  device  that  can  directly  measure  and 
display  the  probability  of  a "l"  on  a signal  line.  The  second  would  be 
a processor  that  can  insert  a number  of  probes  into  a circuit  and  mea- 
sure the  signal  probabilities  of  the  lines  and  process  the  information 
to  relate  to  the  particular  application. 

The  first  monitor  has  been  designed,  and  a prototype  is 
being  exercised.  The  second  monitor  is  being  considered  and  would  be 
implemented  on  an  IBM  System/7  computer.  These  tools  should  be  useful 
augmentations  to  the  other  monitors  described  above  and  should  prove 
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valuable  in  several  applications,  including  test  generation  and  actual 
circuit  testing. 


References 

1.  Burnett,  C.  J.  and  E.  G.  Coffman,  Jr.,  "A  Study  of  Interleaved 
Memory  Systems,"  Proc . AFIPS  Spring  Joint  Computer  Conference,  36, 
AFIPS  Press,  Montvale,  N.J.,  1970,  pp . 467-474. 

2.  Losq,  J.,  "Effects  of  Failures  on  Performance  of  Gracefully  Degrad- 
able Systems,"  Tech.  Note  No.  103.  Digital  Systems  Laboratory, 
Stanford  University,  Stanford,  Calif.,  Dec  1976. 

3.  Losq,  J.,  "Efficiency  of  Compact  Testing  for  Sequential  Circuits," 
Tech.  Note  No.  104,  Digital  Systems  Laboratory,  Stanford  University, 
Stanford,  Calif.,  Dec  1976. 

4.  Parker,  K.  P.  and  E.  J.  McCluskey,  "Analysis  of  Logic  Circuits  with 
Faults  Using  Input  Signal  Probabilities,"  IEEE  Trans,  on  Computers, 
C-24 , 5,  May  1975,  pp.  573-578. 

5.  Schnurmann,  H.  D.,  E.  Lindbloom,  and  R.  G.  Carpenter,  "The  Weighted 
Random  Test  Pattern  Generator,"  IEEE  Trans,  on  Computers,  C-24,  7, 
Jul  1975,  pp.  695-700. 


C.  Project  7181.  THE  2000  TERMINAL  COMPUTING  SYSTEM 

Principal  Investigator:  V.  Cerf 

Staff:  W.  Warren 

1 . Objective 

A study  has  been  undertaken  to  investigate  the  possibilities , 
requirements,  capabil ities , limitations,  and  advantages  of  a computing 
system  serving  a very  large  number  of  users. 


2 . Current  Status  of  Work 

Computer  and  communication  networks  linking  potential  users 
over  a large  geographic  area  make  possible  a computing  system  with  an 
active  user  population  of  several  thousand.  Functional,  rather  than 
general  purpose,  multiprocessing  may  make  such  a system  possible.  It 
could  realize  a number  of  economies  of  scale  and  specialization  and 
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could  provide  reliability  through  modularity  and  redundancy.  We  have 
investigated  the  fundamental  architectural  requirements,  limitations, 
functional  organization,  and  specialization  of  these  systems,  the  soft- 
ware requirements  in  a network  environment,  and  the  functional  charac- 
teristics of  the  user  interface  (both  hardware  and  software). 

In  our  investigation  of  system  architectures  feasible  for  use 
in  a 2000  terminal  computing  system,  we  have  selected  a distributed  net- 
work of  fully  interconnected  microprocessors.  The  microprocessors  will 
be  individually  allocated  to  the  functional  tasks  of  (1)  terminal  I/O 
modules,  (2)  common  memory  modules,  and  (3)  the  general-processing  mod- 
ules. Each  module  will  contain  local  memory,  processing  power,  and  a 
communication  interface  that  will  provide,  via  micro-watt  radio,  a com- 
plete interconnection  network  [1 1 . This  interface  will  handle  module- 
module  {jacket  communications  protocol  based  on  the  carrier  sense  multi- 
ple-access technique  developed  for  the  ALOHA  system  [2,3]. 

One  feature  of  this  system  will  be  the  fast  response  time  seen 
by  the  terminal  users.  Each  terminal  connected  to  the  system  will  have 
its  own  dedicated  terminal-1/0  module  attending  to  line  control,  format- 
ting, and  echoing  [4].  Another  significant  feature  would  be  the  poten- 
tial for  hardware  reliability,  implicit  when  multiple  copies  of  each 
hardware  structure  exist.  Such  a system,  with  appropriately  designed 
software,  would  be  able  to  cope  with  the  complete  failure  of  one  or  more 
hardware  modules,  resulting  in  a minor  degradation  of  response  time  as 
seen  by  the  user.  This  type  of  reliability  will  be  obtained  through  the 
same  methods  and  techniques  used  by  Pluribus  [5] — duplication  of  hard- 
ware, isolation  between  independent  hardware  modules,  and  a software 
system  that  performs  self-consistency  checking  and  recovery. 

We  have  investigated  the  problems  of  recovering  error-free  in- 
formation when  a faulty  hardware  module  is  detected.  If  an  error  is  de- 
tected at  some  point  during  the  processing  of  a task,  it  may  be  necessary 
to  amputate  the  faulty  module,  reconfigure  the  hardware  environment,  and 
rerun  the  task  that  failed.  The  rapid  and  systematic  restoration  of 
service  after  an  error  or  malfunction  would  be  a major  design  goal.  The 
state  of  the  system  (including  system  and  user  data  structures  that  be- 
come invalid)  must  be  restored  to  a previously  consistent  state  so  that 
the  task  can  be  run  again. 
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A system  architecture  was  developed  [6]  to  permit  systematic 
recovery  after  the  occurrence  of  an  error  and  the  restoration  of  erro- 
neous data.  Algorithms  were  developed  to  perform  state  restoration, 
thereby  guaranteeing  that  a computer  system  and  its  associated  data 
bases  and  communication  transactions  will  be  restored  to  an  operational 
state  within  a given  time  and  cost  bound  after  detection  of  a system 
failure.  The  model  used  for  these  algorithms  provides  for  recovery  of 
the  program  state,  revoking  of  erroneous  outputs,  and  restoration  of 
external  data  structures. 

This  research  considered  the  optimization  of  a specific  soft- 
ware strategy — the  rollback  and  recovery  strategy,  within  the  framework 
of  a graph  model  of  program  flow  that  encompasses  communication  inter- 
faces and  data-base  transactions.  Algorithms  were  developed  to  optimize 
the  placement  of  dynamic  recovery  checkpoints . A method  was  presented 
for  statically  precomputing  a set  of  optimal  decision  parameters  for  the 
associated  program  model  and  a run-time  technique  for  dynamically  deter- 
mining the  optimal  placement  of  program-recovery  checkpoints . 
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III.  INTEGRATED  CIRCUITS 


A.  Project  4606.  PRECISION  CHEMICAL  MACHINING  OF  S INGLE -CRYSTAL 

SILICON 

Principal  Investigator:  J.  B.  Angell 

Staff:  P.  Barth 

1 . Objective 

This  project  was  established  to  perfect  techniques  for  fabri- 
cating thin  films  of  single-crystal  silicon  (5  p or  less),  using  chemi- 
cal preferential  etching.  This  work  has  resulted  in  a new  dielectric 
isolation  technique  for  silicon  that  has  several  advantages  over  other 
isolation  techniques.  The  new  method  allows  greater  device  density  than 
does  junction  isolation  in  bulk  silicon  and  maintains  a functional  ad- 
vantage over  silicon-on-sapphire  by  providing  both  MOS  and  bipolar  cap- 
abilities . 

2 . Current  Status  of  Work 

The  work  on  this  project  has  been  very  successful.  Although 
JSEP  support  was  terminated  on  31  August  1976,  the  research  is  continu- 
ing under  new  support  provided  by  a contract  from  the  Naval  Electronics 
Laboratory  Center.  A paper  entitled  "Dielectrically  Isolated  Integrated 
Circuits  Made  with  Preferential  Etching"  was  accepted  for  presentation 
at  the  IEEE  International  Solid-State  Circuits  Conference  on  18  February 
1977;  because  this  paper  summarizes  the  high  points  of  the  work  on  this 
project,  it  is  included  as  the  final  report. 


DIELECTRICALLY  ISOLATED  INTEGRATED  CIRCUITS 
MADE  WITH  PREFERENTIAL  ETCHING 

Phillip  W.  Barth,  James  B.  Angell 

A new  dielectric  isolation  process  based  on  chemical  preferential 
etching  of  silicon  has  resulted  in  reproducible  single-crystal  layers 
less  than  5 pm  (±lpm)  thick  etched  from  a heavily  doped  starting  wafer 
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with  a lightly  doped  epilayer.  The  film  is  separated  into  islands  by 
anisotropic  etching  and  insulated  by  SiO,,,  Si.^Nj,  and  polycrystalline 
silicon  (poly-Si)  to  form  such  structures  as  the  NPN  transistor  in  Fig. 
1.  It  has  the  following  advantageous  features: 


• complete  dielectric  isolation,  eliminating  most  parasitics 

• an  embedded  conducting  layer  of  poly-Si,  which  serves  both 
as  an  additional  layer  of  interconnection  and  as  a lead  to 
the  underside  of  silicon  islands 

• device  geometries  smaller  than  those  achievable  with  junc- 
tion isolation 

• self-alignment  advantages  similar  to  those  of  oxide-wall 
isolation  schemes  (such  as  LOCOS , ISOPLANAR) 

• both  MOS  and  bipolar  capability  for  a functional  advantage 
over  SOS 


Fig.  1.  DIELECTRICALLY  ISOLATED  NPN  TRANSISTOR  FORMED  BY  TRI-POLY 
PROCESS . 
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The  steps  in  the  dielectric  isolation  process,  from  the  starting 


wafer  to  a final  bipolar  transistor  structure,  are  described  below  with 
reference  to  Figs.  2,  3,  and  4. 


(1)  Figure  2:  a 5 pm  thick  epitaxial  layer  is  grown  on  a 

(100\  oriented  silicon  wafer  heavily  doped  with  arsenic 
(3  X lO^/cm^).  A thermal  oxide  is  grown  on  this  epi- 
layer,  and  a layer  of  poly-Si  approximately  200  um  thick 
is  deposited  as  a "handle"  for  later  processing. 

(2)  The  wafer  is  preferentially  etched  in  a mixture  of  1 HF: 
3HN0g : 8CH3COOH , which  removes  the  heavily  doped  starting 
substrate  but  stops  etching  at  the  lightly  doped  epilayer 
[ll  . This  etching  step,  unlike  mechanical  lapping,  is 
tolerant  of  wafer  warpage  commonly  encountered  during 
deposition  of  the  thick  poly-Si  layer  [2] , and  its  stop- 
ping point  is  more  easily  controlled  than  with  mechani- 
cal lapping. 
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Fig.  2.  STAGES  IN  TRI-POLY  PROCESS, 
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(3)  A masked  etch  with  KOH  establishes  islands  of  silicon  which 
will  later  become  active  device  regions. 

(4)  A silicon-nitride  layer  is  deposited  over  the  surface;  this 
layer  will  later  act  as  both  an  insulator  and  an  etch-stop 
region.  Contact  holes  are  defined  in  this  layer,  and  a 1 
pm  thick  poly-Si  layer  is  deposited.  A masking  step  defines 
a pattern  in  this  thin  poly-Si  layer. 

(5)  Figure  3;  an  oxide  layer,  heavily  doped  with  arsenic,  is 
deposited  over  the  thin  poly-Si  layer,  and  then  a thick 
(200  pm)  layer  of  poly-Si  is  deposited  as  a supporting  sub- 
strate for  the  end  structure.  During  this  deposition  step, 
the  As-doped  oxide  acts  as  a diffusion  source  to  lower  the 
resistivity  of  the  1 pm  poly-Si  layer. 

(6)  Figure  4;  the  original  200  pm  layer  of  poly-Si  is  etched 
away,  exposing  the  original  epitaxial  surface.  Three  more 
masking  steps  (the  base  diffusion  need  not  be  a masked  dif- 
fusion) suffice  to  form  the  bipolar  transistor  structure  by 
standard  predeposition  and  drive-in  steps  followed  by  con- 
tact holes  and  aluminum  deposition  and  etching. 


32 


Fig.  4.  TRI-POLY  NPN  TRANSISTOR  WITH  PARTIALLY  SELF-ALIGNED  EMITTER. 

By  returning  to  the  original  epitaxial  surface,  this  tri-poly  pro- 
cess (so  called  because  of  the  three  distinct  poly-Si  layers  present  in 
Fig.  3)  avoids  the  pitted  surface  which  is  present  after  the  preferen- 
tial etching  step.  In  addition,  it  produces  a doping  profile  whose  gra- 
dient is  in  the  right  direction  (doping  is  light  at  the  surface,  heavier 
below) . 

Silicon  nitride  is  used  in  this  process  to  ensure  good  metal  cover- 
age. The  buffered  HF,  used  to  etch  oxide  during  masking  steps,  stops  at 
the  planar  surface  formed  by  Si,  poly-Si,  and  nitride.  This  avoids  a 
"trench"  effect  encountered  during  early  runs  at  the  Si-poly-Si  boundary 
when  only  oxide  was  used  as  an  Insulator. 

Preliminary  work  with  the  tri-poly  process  has  produced  vertical 
NPN  and  lateral  PNP  transistors  which  show  good  minority  carrier  life- 
time and  good  carrier  mobility  in  the  dielectrically  isolated  silicon. 
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and  lateral  PNPs  with  a base- 


NPNs  had  h = 80  with  BV(l|,()  = 22  V , 

width  of  10  .an  on  layout  had  )i  „ =5, 

fe 

A difference  in  structure  is  evident  between  the  transistors  in 
Figs.  1 and  4;  in  botli , the  base  diffusion  can  be  an  unmasked  self- 
aligned  diffusion.  The  transistor  in  Fig.  1 requires  complete  alignment 
of  the  emitter,  and  the  emitter  in  Fig.  4 is  partially  aligned  by  using 
the  nitride  wall  as  one  boundary.  Figure  4 also  contains  a parasitic 
MOSFET  (emitter  = source,  base  = channel,  collector  = drain  and  gate) 
which  leads  to  emitter-collector  leakage.  Efforts  are  being  made  to 
eliminate  this  parasitic  MOSFET  action  to  obtain  the  advantages  of  the 
partially  self-aligned  emitter. 

Figure  5 illustrates  the  good  layout  density  available  with  the 
tri-poly  process.  Pictured  is  an  inverter  stage  of  a complementary 


Fig.  5.  SECTION  OF  COMPLEMENTARY  TRANSISTOR  LOGIC 
RING  OSCILLATOR . Lateral  PNP:  base  = Vn , emit- 

ter = i collector  = V^n.  Vertical  NPN : base  = 
V , emitter  = ground,  collector  = V 
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transistor  logic  ring  oscillator,  similar  to  an  I L gate  with  the  fol- 
lowing variations  on  the  I^L  structure:  (a)  the  PNP  base  lead  goes  to 

a separate  voltage  line,  (b)  the  NPN  action  is  downward  rather  than  up- 
ward, and  (c)  the  NPN  is  Schottky  clamped.  This  gate  requires  two  extra 

2 

aluminum  lines  (V  and  ground)  as  compared  to  an  I L gate  in  bulk  sil- 
B 

icon,  and  extra  space  is  taken  up  by  the  Schottky  diode;  however,  the 

2 

area  is  only  1.6  times  that  of  an  I L gate  using  the  same  layout  rules. 

2 

A straight  reproduction  of  an  I L inverter  using  the  tri-poly  process 

(including  an  extra  line  for  ground)  would  have  an  area  0.6  times  that 
2 

of  a bulk  I L inverter. 
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B.  Project  5012.  IC  PROCESS  DESIGN  AND  COMPUTER  MULTILAYER 

SEMICONDUCTOR  DEVICE  ANALYSIS 

Principal  Investigator:  R.  W.  Dutton 

Staff:  D.  B.  Estreich 


1 . Ob jectlves 


Computer  methods  were  developed  to  predict  impurity  profiles 
after  processing,  and  physical  parameters  related  to  impurity  diffusion 
were  calculated. 

Improved  models  will  be  developed  for  the  prediction  and  sim- 
ulation of  multilayer  semiconductor  device  behavior  based  on  layout  and 
device-processing  history. 

2.  Current  Status  of  Work 

During  this  past  reporting  period,  the  following  major  accom- 
plishments have  been  realized: 
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• Improved  parameter  selection  for  a new  user-oriented 
integrated  injection  logic  (I^L)  macromodel  that  models 
current  gain  and  injector  transport  efficiency  falloff 
effects 

• evaluation  and  design  of  second  generation  of  special 
complementary  MOS  (CMOS)  latch-up  path  test  structures 
to  allow  separation  of  the  relevant  parameters  associ- 
ated with  PNPN  latch  paths 

The  I~L  macromodel  and  CMOS  latch-up  work  are  a continuation  of  the  re- 
search described  in  the  last  Status  Report. 

2 

Previous  work  on  modeling  1 L has  been  reported  [1,2];  how- 
ever, the  models  presented  were  applicable  only  to  low-current  operation. 

2 

Recent  emphasis  has  focused  on  modeling  those  aspects  of  I L that  pre- 
sent severe  limitations  to  I~L  performance  and  operational  bounds  [3,41. 
This  requires  the  modeling  of  both  NPN  current  gain  falloff  resulting 
from  base  resistance  (and  emitter  resistance  in  certain  instances)  and 
injector  PNP  current  transport  degradation  at  low-  and  high-current  op- 
erating levels.  Parameter  measurement  schemes  have  been  developed  and 
applied  for  base  resistance  and  upward  transit-time  determination,  in- 
cluding proper  partitioning  of  multiple  fan-out  gates.  With  the  inclu- 
sion of  these  effects  in  the  macromodel,  excellent  agreement  can  be 

achieved  between  computer  simulation  and  experimental  results  for  all 
2 

I L operating  bias  levels  . 

Recent  work  has  also  involved  the  study  of  latch-up  mechanisms 
in  CMOS  integrated  circuits  and  related  merged  structures  . Two-dimen- 
sional analysis  is  required  for  adequately  modeling  latch-up  action  in 
integrated  circuits.  It  has  been  demonstrated  that  accurate  accounting 
for  lateral  spreading  resistance  is  crucial  in  predicting  latch-up  con- 
ditions. In  addition,  field-aided  lateral  PNP  current  gain  enhancement 
is  being  studied.  New  test  structures  have  been  designed  to  examine 
these  factors . 
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IV.  SOLID  STATE 


A.  Project  5111.  TRANSPORT  PROPERTIES  OF  A^Ga^^s,  A1  Ga  Sb , AND 

In  Ga,  As  SINGLE  CRYSTALS  X 1_X 

x 1-x 

Principal  Investigator:  G.  L.  Pearson 

Staff:  K.  Jew 

1 . Objective 

The  purpose  of  this  project  is  a systematic  study  of  the  trans- 
port of  free  carriers  in  these  ternary  systems . The  electronic  parameters 
are  free  carrier  mobility,  resistivity,  minority  carrier  diffusion  length  , 
activation  energy,  and  trap  energy  levels.  The  distribution  coefficients 
of  selected  doping  impurities  as  a function  of  x are  significant  growth 
parameters.  The  results  obtained  should  prove  useful  in  the  design  oi 
photoluminescent  devices,  injection  lasers,  solar  cells,  infrared  detec- 
tors, ar.d  bulk-effect  microwave  devices. 


2 . Current  Status  of  Work 

Epitaxial  layers  of  Al^Ga^  ^As  were  grown  by  liquid  phase  epi- 
taxy techniques  on  [100]  GaAs  substrates.  The  melt  solutions  were  pre- 
pared with  6N  gallium,  polycrystalline  GaAs,  and  aluminum  metal.  Dopant 
impurities  were  Sn  for  n-type  layers  and  Ge  for  p-type  layers.  To  reduce 
unintentional  impurities,  the  gallium  melt  was  initially  baked  for  more 
than  12  hours  at  800 °C  in  flowing  high-purity  hydrogen.  After  cooling 
to  room  temperature,  substrates  were  inserted  and  source  materials  were 
added.  Growth  was  carried  out  following  saturation  at  800°C  through  1 °C 
of  supercooling  before  bringing  the  substrate  into  contact  with  the  melt. 
The  cooling  rate  was  0.2°C/min  and  growth  terminated  at  780 °C. 

During  this  investigation,  two  crystal-growth  procedures  have 
been  used.  In  the  first,  the  same  melt  was  reused  for  several  growth 
runs.  Because  the  aluminum  content  in  successive  grown  layers  decreases 
as  a result  of  depletion  of  A1  from  the  melt,  a range  of  ternary  compo- 
sitions can  be  obtained.  In  the  second  method,  a fresh  melt  was  prepared 
before  each  growth  run.  Most  of  the  epitaxial  layers  were  grown  from 
reused  melts.  This  method  was  discontinued,  however,  when  it  was  found 
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that  (1)  precise  growth  data  cannot  be  obtained  because,  in  a series  of 
growths,  materials  are  depleted  from  the  melt  and  the  exact  melt  compo- 
sition is  unknown  after  the  initial  growth,  and  (2)  the  reused  melts 
become  contaminated  as  a result  of  oxidation  of  the  A1  during  repeated 
cyclings  of  the  system  between  room  and  growth  temperatures,  making  it 
increasingly  difficult  to  grov'  good  quality  crystals.  Currently, 
Al^Ga^^As  layers  are  grown  froi,-’  fresh  melts  only.  These  samples  are 
expected  to  yield  more  reliable  results  for  dopant  and  A1  distribution 
coefficients  and  for  electrical  properties. 

Photoluminescence  measurements  were  made  at  77°  and  300 °K  on 
samples  of  (nlAl^Ga^  ^AsiSn  with  composition  0.0<x<0.50.  The  band- 
gaps  and  alloy  compos itions  were  determined  from  the  photoluminescence 
peak . 

Electrical  measurements  of  (nlAl^Ga^  ^AsiSn  samples  were  made 
using  Van  der  Pauw  techniques.  Ohmic  contacts  were  formed  by  Au-Ge  eu- 
tectic and  alloying  in  hydrogen  at  450 °C.  Carrier  concentration,  elec- 
tron mobility,  and  resistivity  were  measured  from  77°  to  300°K.  From 
the  temperature  dependence  of  the  carrier  concentration,  the  thermal 
activation  energy  of  the  Sn  donor  has  been  determined  for  Al  Ga  As  in 

the  0.0  < x < 0.50  range.  Results  depend  on  the  method  used  to  grow 
the  epitaxial  layers.  For  (nlAl^Ga^  ^As :Sn  crystals  of  similar  compo- 
sitions and  carrier  concentrations,  the  electron  mobility  of  crystals 
grown  from  fresh  melts  is  generally  higher  when  compared  to  crystals 
grown  from  reused  melts  as  a result  of  the  decrease  in  crystalline  qual- 
ity. 

Current  study  is  directed  toward  LPE  (plAl^Ga^  ^As:Ge  layers 
with  compositions  0 < x < 1.0.  Van  der  Pauw,  photoluminescence,  and 
scanning  electron  microscope  techniques  are  employed. 
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B.  Project  5244.  STUDIES  OF  HIGH-TRANSITION  SUPERCONDUCTORS  SUCH  AS 

V3Si  AND  Nb3Sn  AND  DEPTH  PROFILING  OF  A1  Ga  As- 
GaAs  HETERO, JUNCTIONS  X ~X 

Principal  Investigator:  W.  E.  Spicer 

Staff:  I.  Lindau , J.  Miller,  C.  M.  Garner 

1 . Objective 

This  project  is  involved  in  two  separate  areas  of  research. 

The  first  is  to  explore  the  electronic  structures  of  type-II  (hard)  su- 
perconducting alloys,  particularly  the  structure  of  Nb^Sn , to  determine 
such  parameters  as  valence  bandwidth  which  is  important  to  the  theoret- 
ical models  of  superconductors.  Emphasis  is  directed  toward  the  prepa- 
ration of  atomically  clean  Nb  surfaces,  followed  by  evaporation  of  Sn 
as  an  adlayer  and  annealing  to  form  Nb^Sn. 

The  second  is  to  study  the  chemical  interface  between 

Al^Ga^  ^As  and  GaAs  in  this  hetero junction  system,  particularly  the 

chemical  width  of  the  junction  and  surface  morphology  of  the  interface 

for  various  aluminum  concentrations  (0.3  < X < 0.85).  The  Schottky 

barrier  Au-Af  Ga,  As  interface  is  also  to  be  examined, 
x 1 — x 

2 . Current  Status  of  Work 

a.  High-Transition  Superconductors  (I.  Lindau,  J.  Miller) 

Earlier  work  on  a Nb  single  crystal  indicated  that  the 
primary  contaminants  of  the  Nb  surface  were  carbon,  nitrogen,  and  oxygen. 
As  a result,  several  Nb  foils  were  purchased  to  investigate  new  methods 
for  cleaning  the  samples.  To  determine  the  origin  of  the  surface  con- 
tamination, we  elected  to  depth  profile  the  carbon  and  oxygen  signals. 

The  Nb  foils  of  better  than  99.8  percent  bulk  purity  were  subjected  to 
continuous  cleaning  by  argon-ion  sputtering.  The  Auger  signals  of  oxy- 
gen and  carbon  were  constantly  monitored  during  the  sputtering,  and  a 
clean  surface  was  established  with  the  carbon  and  oxygen  signals  falling 
below  the  noise  level.  The  oxygen  contamination  reappeared  immediately 
after  sputtering,  however,  or  when  the  argon  beam  energy  was  reduced  be- 
low 500  eV . 

Measurements  of  the  partial  pressure  of  0^  and  CO  with  a 
mass  spectrometer  indicated  that  all  of  the  contamination  could  not  stem 
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from  the  background  gases  present  in  the  chamber.  One  mechanism  proposed 
to  explain  the  almost  immediate  reappearance  of  the  contaminants  was  dif- 
fusion from  the  bulk. 

In  other  work,  successful  removal  of  carbon  was  accomplished 
by  heating  the  sample  to  1700°C  for  three  hours.  Light  sputtering  removed 
the  nitrogen,  but  most  of  the  oxygen  signal  remained. 

b.  Depth  Profiling  (C.  M.  Garner) 

Auger  profiles  have  been  made  on  Af^Ga^  ^As-GaAs  hetero- 
junctions [grown  by  abrupt  liquid  phase  epitaxy  (LPE)  techniques]  for 
x = 0.3,  0.43,  0.60,  and  0.85.  The  interface  width  (90  to  10  percent  A S 
Auger  signal)  was  found  to  be  ~100  \ for  samples  grown  at  750°C  and  130  A 
for  samples  grown  at  800°C.  The  growth  mechanisms  and  dependence  of  the 
interface  width  on  growth  temperature  and  time  are  being  investigated. 

Additional  studies  have  included  the  InGaAsP-InP  interface 
of  samples  grown  by  LPE  and  VPE . The  LPE  samples  were  found  to  behave  as 
predicted  by  the  growth  conditions;  however,  the  VPE  InP  layer  displayed 
an  unusual  composition  with  depth  (it  appeared  that  the  Ga  had  diffused 
out  of  the  InGaAsP  layer  into  the  InP).  Further  work  will  be  directed 
toward  understanding  this  effect. 

Auger  profiles  of  the  Au-Af^Ga^  ^As  interface  revealed  the 
presence  of  an  oxide  layer  in  samples  that  had  been  etched  prior  to  being 
placed  in  the  Au  evaporation  chamber;  this  layer  is  approximately  the 
same  thickness  as  the  native  oxide  layer.  Current  work  is  directed  at 
removing  the  oxide  layer  and  evaporating  Au  in  situ  to  avoid  the  forma- 
tion of  an  oxide  layer.  The  electrical  properties  of  these  contacts  will 
then  be  studied  to  understand  the  effect  of  the  oxide. 
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V. 


RADIOSCIENCE 


A.  Prolect  3172.  DIGITAL  IMAGE  PROCESSING 

Principal  Investigator:  R.  N.  Bracewell 

Staff:  J.  Mandeville 

1 . Objective 

The  purpose  of  this  project  is  to  evaluate  and  improve  algo- 
rithms for  digital  image  processing  with  a priori  information. 

2.  Current  Status  of  Work 

This  project  will  become  active  on  1 January  1977. 


B.  Project  3476.  DEVELOPMENT  OF  AN  AUDIO-TRACKING  FILTER/ANALYZER 

Principal  Investigator:  R.  A.  Helliwell 

Staff:  E.  Paschal 

1 . Objective 

Variable  frequency  narrowband  VLF  signals  in  the  magnetosphere 
are  difficult  to  analyze.  At  present,  we  record  VLF  signals  at  field 
stations  and  process  the  recordings  at  Stanford,  using  a spectrum  ana- 
lyzer to  generate  photographic  spectral  records.  It  is  difficult,  how- 
ever, to  determine  the  amplitude  of  individual  signals  in  the  filmed 
spectra.  One  solution  is  to  use  a frequency-tracking  filter  that  can 
follow  the  changing  frequency  of  an  individual  signal  and  separate  it 
from  other  signals  present  in  the  tape-recorded  data  so  that  its  ampli- 
tude as  a function  of  time  can  be  measured.  The  purpose  of  this  project 
is  to  develop  such  a filter,  capable  of  tracking  signals  over  the  range 
from  50  Hz  to  30  kHz.  The  filter  will  have  selectable  bandwidths  from 
10  Hz  to  1 kHz.  wide  and  will  be  capable  of  tracking  signals  that  change 
frequency  at  rates  up  to  20  kHz/sec. 
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Current  Status  of  Work 


The  design  of  the  system  is  approximately  90  percent  complete, 
and  the  construction  is  70  percent  finished.  We  estimate  that  another 
two  man-months  will  be  required  to  complete  the  project.  Because  the 
project  leader  (E.  Paschal)  will  be  working  in  the  Antarctic  from  Octo- 
ber 1976  until  March  1977,  the  project  should  be  finished  sometime  in 
May  1977. 

The  tracking  filter  uses  a local  oscillator  and  mixer  to  trans- 
late the  input-signal  spectrum  into  the  passband  of  an  IF  filter  with  a 
center  frequency  of  100  kHz.  The  output  of  the  filter  is  then  fed  to  a 
frequency  discriminator  whose  output  controls  the  LO  frequency  to  cause 
the  system  to  follow  the  input-signal  component  as  it  changes  frequency. 
The  tracking  filter  system  can  be  divided  into  two  parts — the  IF  filter 
and  discriminator  assembly  and  the  local  oscillator  and  track-control 
assembly . 


a . IF  Filter  and  Discriminator  Assembly 

The  IF  filter  uses  conventional  analog  techniques  to  fil- 
ter the  input  signal;  actually,  several  filter  elements  are  used  in  se- 
ries. The  input  signal  at  100  kHz  is  run  through  a three-pole  Butterworth 
filter  with  a bandwidth  of  3 kHz  and  is  then  mixed  with  a 70  kHz  tone. 
The  30  kHz  difference  spectrum  is  then  filtered  with  another  three-pole 
filter  with  a 1 kHz.  bandwidth,  and  the  output  is  translated  down  to  10 
kHz.  This  process  is  repeated  with  successive  filters  with  bandwidths 
of  300,  100,  and  30  Hz  and  finally  with  a 10  Hz  bandwidth  filter  centered 
at  300  Hz.  System  bandwidths  of  10,  30,  100,  300  Hz,  and  1 kHz  are  ob- 
tained by  selecting  the  output  of  the  appropriate  filter.  Each  filter 
has  a matching  discriminator,  and  the  output  of  the  selected  discrimina- 
tor is  used  to  control  the  local  oscillator.  The  IF  filters  and  discrim- 
inators employ  temperature-compensated  L-C  resonators. 

This  multiple-conversion  approach  makes  it  possible  to 
build  a variable  bandwidth  filter  with  a given  (input)  frequency  without 
resorting  to  crystal  resonators  which  would  be  required  if  all  the  fil- 
tering were  to  be  done  at  100  kHz.  The  spurious  response  to  out-of-band 
signals  should  also  be  less  than  with  a crystal  filter.  The  IF  filter 
and  discriminator  assembly  are  nearly  one-half  completed. 
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Local  Oscillator  and  Track  Control 


This  part  of  the  system  uses  all  digital  circuits  and  has 
been  the  most  interesting  segment  of  the  design.  To  track  a signal,  the 
local  oscillator  is  set  to  such  a frequency  that  the  signal,  when  it  ap- 
pears, will  be  translated  into  the  passband  of  the  IF  filter.  The  LO 
frequency  runs  from  100  to  130  kHz.  It  is  necessary  to  set  the  initial 
frequency  of  the  LO  with  sufficient  accuracy  so  that  the  translated  sig- 
nal is  within  the  passband  of  the  IF  filter;  with  a filter  bandwidth  of 
10  Hz,  this  means  that  the  oscillator  must  be  set  within  ±5  Hz  of  the 
correct  frequency.  Because  such  accuracy  at  a frequency  of  100  kHz  was 
deemed  beyond  the  ability  of  an  analog  oscillator,  a digital  frequency 
synthesizer  has  been  designed  for  the  LO. 

The  synthesizer  uses  a novel  form  of  direct  digital  syn- 
thesis. Its  output  is  a square  wave  at  a frequency  from  100  to  130  kHz. 
The  litter  in  the  square  wave  is  approximately  3 nsec  rms  and,  as  a re- 
sult, spurious  synthesizer  signals  are  at  least  60  dB  below  the  funda- 
mental. (A  spurious  LO  signal  appears  in  the  system  as  a spurious  filter 
response  and  must  be  kept  as  small  as  possible.)  The  synthesizer  can  be 
controlled  in  1 Hz  steps  and  has  an  accuracy  of  0.1  Hz.  CMOS  logic  is 
used,  and  the  synthesizer  consumes  approximately  5W  of  power.  Its  design 
is  apparently  new  and  may  be  patented . 

Using  a digital  frequency  synthesizer  for  the  local  os- 
cillator also  simplifies  the  frequency  and  tracking-control  circuits. 
The  initial  frequency-set  and  the  frequency  and  slew-rate-limit  circuits 
are  all  digital.  The  construction  of  the  local  oscillator  and  track- 
control  circuits  is  approximately  95  percent  complete. 


C.  Projects  3606  and  4504.  MEASUREMENTS  OF  METEOROLOGICAL  PARAMETERS 

IN  THE  LOWER  TROPOSPHERE 

Principal  Investigators:  A.  M.  Peterson,  V.  R.  Eshleman, 

A . T . Waterman , Jr  . , M . Frankel 
Staff:  N.  Bhatnagar,  S . J . Wernecke , R.W.Lee,R.D. 

Fleming 

1 . Ob Jective 

The  objective  of  this  research  is  the  investigation  of  atmo- 
spheric parameters  and  phenomena  in  the  lower  troposphere  based  on  two 
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different  but  complementary  techniques — radio  acoustic  sounding  (HASS) 
and  radio  propagation  in  the  S-  and  K-bands.  The  goals  of  each  method 
are  as  fo 1 lows . 

a.  HASS  (Project  3606) 

• Measuring  real-time  vertical  temperature  and  horizontal 
wind  profiles  over  a 1 km  (or  greater)  range,  using  both 
RASS  and  conventional  techniques. 

• Obtaining  temperature  and  horizontal  wind  profiles  over 
extended  periods  of  time  and  through  varying  meteorolog- 
ical conditions. 

• Monitoring  atmospheric  phenomena  continuously  (with  sam- 
pling periods  of  less  than  5 minutes). 

• Studying  the  interaction  of  electromagnetic  and  acoustic 
waves  in  a stochastic  atmosphere  characterized  by  turbu- 
lence, mean  winds,  and  temperature  gradients. 

b.  Tropospheric  Radio  Propagation  (Project  4504) 

• Measurement  of  wind  velocity  in  the  common  volume  of  a 
transhorizon  troposcatter  propagation  path  by  means  of 
doppler  techniques. 

• Probing  the  structure  of  winds  and  turbulence  in  the 
lower  troposphere  under  various  conditions  of  atmo- 
spheric stability  and  instability. 

2 . Current  Status  of  Work 

a.  RASS  (Project  3606) 

In  the  Stanford  Radio  Acoustic  Sounding  System  (RASS), 
an  electromagnetic  signal  is  made  to  scatter  from  a moving  acoustic 
pulse  train.  Under  a Bragg-sca tter  condition,  maximum  electromagnetic 
scattering  occurs.  The  scattered  radio  signal  contains  the  signature 
of  the  atmosphere  as  a function  of  the  acoustic-pulse  position. 

In  the  analysis  of  the  RASS  to  date,  such  atmospheric 
conditions  as  turbulence,  humidity,  mean  temperature,  and  mean  wind 
fields  have  been  ignored.  In  this  study,  system  performance  is  assessed 
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in  a real  atmosphere.  The  only  assumption  marie  is  that  the  acoustic 
wave  train  (not  the  EM  wave)  is  affected  by  stochastic  perturbations  in 
the  atmosphere. 

Coherency  of  vertical  acoustic-wave  propagation  is  de- 
scribed through  a perturbation-theoretic  method  and  Feynman's  diagram- 
matic technique.  One  of  the  most  important  attributes  of  this  a-  'lysis 
is  that  it  systematically  and  explicitly  accounts  for  multiple  scatter- 
ing of  acoustic  waves  in  the  presence  of  atmospheric  fluctuations.  The 
coherency  results  are  then  applied  to  evaluate  the  strength  of  the  scat- 
tered electromagnetic  signal  from  the  acoustic  pulse  train,  while  taking 
into  account  the  presence  of  turbulence,  mean  temperature,  and  mean  wind 
f ields  . 

From  this  study,  we  conclude  that,  for  acoustic  pulses 

with  carrier  frequencies  below  a few  kilohertz  propagating  under  typical 

atmospheric  conditions,  turbulence  has  little  effect  on  the  strength  of 

the  received  radio  signal  at  heights  up  to  a few  kilometers.  This  result 

implies  that  focusing  of  RF  energy  by  the  acoustic  wavefronts  is  primar- 

-2 

ily  a function  of  sound  intensity  which  decreases  as  x , where  x is 
the  altitude. 

The  effect  of  mean  vertical  wind  and  mean  temperature  on 
the  strength  of  the  received  signal  is  also  shown  to  be  insignificant. 
Mean  horizontal  winds,  however,  shift  the  focus  of  the  reflected  elec- 
tromagnetic energy  from  its  origin,  resulting  in  a decrease  in  received 
signal  level  when  a monostatic  RF  system  is  used.  In  a bistatic  radar 
configuration  with  space-diversified  receiving  antennas,  the  shifting  of 
the  acoustic  pulse  makes  possible  the  remote  measurement  of  the  horizon- 
tal wind  component. 

These  theoretical  analyses  have  been  supported  experimen- 
tally. A high-frequency  RASS  operating  at  an  acoustic  frequency  of  1 kHz 
and  a radio  frequency  of  440  MHz  was  constructed  at  the  Stanford  Research 
Institute.  Consistent  with  theoretical  predictions,  this  system  was 
able  to  provide  both  doppler  data  and  horizontal  wind  profiles  over  a 
range  of  1 km. 

These  experimental  results  came  from  a brief  endeavor  to 
support  our  theoretical  studies.  Future  experiments  could  extend  the 
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profiling  range  and  versatility  of  the  high-frequency  RASS  ; it  could 
become  a valuable  tool  for  the  remote  real-time  measurement  of  atmo- 
spheric temperature,  winds,  and  wind-shear.  Through  further  research, 
we  hope  to  provide  a mobile  system  capable  of  producing  these  data  at 
any  location  desired  by  a user. 

During  the  next  reporting  period,  we  intend  to  complete 
our  theoretical  research  and  submit  it  for  publication. 

b.  Tropospheric  Radio  Propagation  (Project  4504) 

As  outlined  in  the  last  status  report,  the  wind -measur- 
ing technique  is  to  involve  the  simultaneous  reception  of  3-GHz  CW  sig- 
nals on  each  of  two  narrow  beams  from  a transmitter  located  100  miles 
away  (beyond  1 ine-of-s ight ) . The  two  beams  are  aimed  on  either  side  of 
the  great-circle  bearing  toward  the  transmitter  by  small  amounts  (1°  or 
less).  If  there  is  a wind  component  transverse  to  the  path  in  the  re- 
gion of  the  common  volume,  it  is  anticipated  that  the  signal  received 
on  one  beam  will  be  doppler  shifted  to  higher  frequencies  and  the  sig- 
nal on  the  other  to  lower  frequencies.  The  magnitude  of  the  doppler 
difference  should  be  proportional  to  transverse-wind  velocity.  The  two 
beams  will  be  formed  during  the  data -reduct ion  procedures  because  the 
receiving  antenna  is  a sampling  array  that  measures  and  stores  on  dig- 
ital tape  the  amplitudes  and  phases  of  the  signal  as  received  on  12  in- 
dividual antennas.  To  obtain  high  resolution  in  azimuth,  the  antennas 
are  arranged  in  a horizontal  linear  array. 

The  linear  array  to  be  used  for  this  purpose  is  a series 
of  12  4-ft  parabolic  dishes  mounted  on  a 70-ft  tower.  The  tower  had 
been  used  in  previous  experiments  in  a vertical  position.  For  the  cur- 
rent work,  it  was  necessary  to  move  the  tower  to  a horizontal  position, 
and  this  move  has  been  accomplished.  The  balance  of  the  equipment  has 
been  checked  and  is  in  operating  condition;  however,  the  measurement 
program  has  been  delayed  because  one  of  the  key  personnel  is  on  leave 
of  absence. 

In  the  meantime,  some  work  has  been  directed  toward  the 
second  of  the  two  objectives,  involving  1 ine-of-s ight  rather  than  trans- 
horizon paths. 
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I).  Project  4214.  INVESTIGATION  OF  UNDERSEA  COMMUNICATION  WITH  ULF 

ELECTROMAGNETIC  WAVES 

Principal  Investigator:  O.  G.  Villard,  Jr. 

Staff:  A.  C.  Fraser-Smith , D.  M.  Bubenik,  A.  Ho 

1 . Objective 

The  feasibility  of  using  ULF  electromagnetic  waves  (frequen- 
cies of  less  than  5 Hz ) from  electric  and  magnetic  dipole  sources  will 
be  investigated  for  undersea/undersea  communication  in  the  approximate 
range  of  0 to  20  km. 

2 . Current  Status  of  Work 

This  is  a new  project;  work  commenced  in  September  1976.  The 
theoretical  expressions  for  the  ULF  electric  and  magnetic  fields  pro- 
duced in  sea  water  by  submerged  vertically  directed  magnetic  and  elec- 
tric dipoles  (with  a sea  floor  present)  have  been  programmed  for  numer- 
ical solution.  Numerical  solution  is  in  progress,  and  exact  data 
concerning  ULF  fields  will  become  available  shortly.  We  will  then  be 
in  a position  to  evaluate  the  feasibility  of  using  these  fields  for 
communication  beneath  the  sea. 
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VI . PLASMA  PHYSICS  AND  QUANTUM  ELECTRONICS 

A.  Project  1328.  VERY  LONG  DELAYED  RADIO  ECHOES 

Principal  Investigator:  F.  W.  Crawford 

Staff:  R.  J.  Vidmar,  T.  L.  Savarino 

1 . Objective 

There  is  evidence  indicating  that  radio  signals  can  occasion- 
ally be  returned  from  the  ionosphere  with  delays  of  the  order  of  tens  of 
seconds,  rather  than  the  few  milliseconds  corresponding  to  direct  re- 
flection. The  objective  of  this  project  is  to  develop  a theoretical 
explanation  of  this  phenomenon. 

2 . Current  Status  of  Work 

Earlier  theoretical  studies  [1]  have  suggested  the  following 
plausible  mechanism  for  long-delayed  echoes.  When  an  ordinary  wave  from 
the  transmitter  reaches  its  reflection  point,  near  the  peak  of  the  F- 
layer,  it  is  partially  converted  into  a longitudinal  plasma  wave.  This 
wave  is  amplified  by  energetic  electrons  precipitating  down  Earth's  mag- 
netic field  lines  and  propagates  with  very  low  group  velocity  («1  km/ 
sec)  over  a distance  of  several  kilometers.  The  longitudinal  plasma 
wave  is  partially  converted  into  a transverse  electromagnetic  wave  as 
it  propagates  down  the  ionospheric  electron-density  gradient,  and  this 
wave  propagates  to  the  receiver. 

Work  was  performed  on  this  project  only  during  the  first  two 
months  of  the  reporting  period.  Effective  1 September  1976,  the  project 
was  transferred  to  other  sponsorship.  The  work  performed  before  termi- 
nation was  directed  toward  an  understanding  of  the  process  by  which  the 
amplified  longitudinal  plasma  wave  is  coupled  to  a transverse  electro- 
magnetic wave  which  can  propagate  out  of  the  ionosphere.  Specifically, 
the  coupling  into  an  ordinary  wave  has  been  studied  in  the  region  w « 

u , near  the  ordinary  wave  reflection  point.  Coupling  between  these 
P 

wave  modes  is  assumed  to  occur  as  a result  of  Earth's  magnetic  field  and 
an  inhomogeneity  in  the  ambient  electron  density.  Even  with  a simple 
model  of  this  amplified  longitudinal  wave,  significant  coupling  into  the 
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ordinary  wave  is  predicted,  confirming  preliminary  estimates  made  in  the 
last  reporting  period.  The  more  complicated  step,  to  develop  a model  of 
the  mechanism  that  couples  energy  into  longitudinal  waves  from  the  elec- 
tromagnetic waves  produced  by  the  ground  transmitter,  is  being  pursued 
under  NSF  sponsorship. 

Reference 

1.  D.  M.  Sears,  "Long  Delayed  Echoes,"  Stanford  University  Institute 
for  Plasma  Research  Report  No.  584,  Stanford  University,  Stanford, 
Calif.,  Nov  1974. 


B.  Hansen  Laboratories.  TWO-PHOTON  RESONANTLY  PUMPED  IR  UP-CONVERTERS 

Principal  Investigators:  S.  E.  Harris,  J.  F.  Young 

Staff:  J.  H.  Newton 

1.  Objective 

The  goals  of  this  project  are  the  development  and  extension  of 
efficient  IR  up-conversion  techniques  in  metal  vapors,  particularly  im- 
age up-conversion. 

2.  Current  Status  of  Work 

During  the  previous  reporting  period , we  demonstrated  the  up- 

conversion  of  2.9  urn  IR  images  to  4550  A in  Cs  vapor.  A power  efficiency 

of  20  percent  with  over  1000  resolvable  spots  was  achieved  using  a pump 

2 

power  of  only  8 kW  focused  in  an  area  of  0.1  cm  . The  Cs  cell  had  an 

17 

active  length  of  0.2  cm  and  was  operated  at  a density  of  1.4  X 10  /cc . 

The  pumping  laser,  Nd: lanthanum  berylate,  has  a natural  two-photon  coin- 

2 2 

cidence  with  the  Cs  6s  S-7s  S transition,  resulting  in  a simple  practi- 
cal system  with  a number  of  potential  applications. 

We  have  since  designed  and  constructed  a new  Cs  IR  imaging 
system  that  should  have  significantly  increased  resolution  and  effi- 
ciency. As  seen  in  Table  1,  this  increase  will  be  accomplished  through 
higher  number  density  Cs  cells  of  larger  area.  For  comparison,  Table  1 
also  lists  the  parameters  for  a LiNbO.^  up-converter  having  the  same  per- 
formance; the  crystal  system  requires  1250  times  more  pump  power. 
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Table  1 


DESIGN  PARAMETERS  FOR  A Cs  AND  LiNbC>3  2.9  p UP-CONVERTER 


Parameter 

Cs 

LiNb03 

Pump  wavelength  (p) 

1.0790 

1 .06 

IR  wavelength  (p) 

2.94 

2 .94 

Sum  wavelength  (A) 

4560 

7800 

Pump  bandwidth  (cm  *) 

0.1 

=4.0 

Pump  pulselength  (nsec) 

50 

50 

Photon  efficiency  ($) 

27.5 

27.5 

Resolution 

io5 

io5 

Pump  power  (W) 

1.44  X 105 

9 X 107 

Cs  density  (cm  3) 

17 

4.44  X 10 

— 

2 

Area  (cm  ) 

3.74 

2 

Length  (cm) 

0.2 

0.46 

As  a result  of  attacks  by  the  Cs  vapor,  our  earlier  cells  have 
had  short  lifetimes  (approximately  24  hours).  The  new  design  Incorpor- 
ates features  that  should  make  it  more  Impervious  to  attack  and  substan- 
tially Increase  the  lifetime. 

The  cell  used  In  the  demonstration  of  the  IR  up-conversion  con- 
sisted of  two  metalized  sapphire  windows  brazed  to  kovar  rings.  The  rings 
were  heliarced  to  Varian  conflat  flanges  that  were  then  bolted  together 
by  means  of  standard  copper  gaskets.  A nickel  side  arm  was  heliarced  to 
the  flanges.  The  spacing  between  the  windows  was  0.2  cm.  The  apparent 
failure  of  the  cell  occurred  at  the  braze  between  the  sapphire  and  kovar 
ring.  As  a result,  we  have  decided  on  a radically  different  cell  design 
(see  Fig.  6)  consisting  of  two  2-inch  in  diameter  sapphire  windows  pressed 
against  an  optically  polished  stainless-steel  plate  0.2  cm  thick.  In  the 
center  of  the  plate  will  be  a hole  1.25-inch  in  diameter  to  produce  a 3/8- 
lnch  seal  around  each  window.  Each  sapphire  window  will  be  pressed 
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Fig.  6.  NEW  Cs  CELL  DESIGN. 

against  the  polished  plate  by  means  of  an  annulus  spring  loaded  by  a stack 
of  Belleville  spring  washers  held  by  the  port  window  flanges  of  the  cylin- 
drical outer  cell.  Two-inch  infrasil  quartz  windows  will  be  rubber  o-ring 
sealed  to  the  port  window  flanges.  The  stainless-steel  plate  will  be  elec 
tron-beam  welded  to  a hollow  cylindrical  stainless-steel  support,  and  the 
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support  will  be  metal  o-ring  sealed  to  the  base  plate.  The  base  plate 
has  a thin-walled  stainless-steel  side  arm  heliarced  to  it . A 0 . 1 cm 
hole  connects  the  hole  in  the  center  of  the  polished  stainless  plate  and 
the  side  arm,  and  a capillary  tube  connects  the  side  arm  to  a hole  in 
the  base  plate.  To  load  the  cell,  a clean  glass  ampoule  of  Cs  is  in- 
serted in  the  side  arm.  Evacuating  the  outer  cell  evacuates  the  space 
between  the  sapphire  windows  and  the  side  arm  via  the  capillary  tube. 

The  ampoule  is  then  crushed  by  flexing  the  thin-walled  side  arm.  When 
the  side  arm  is  heated  , Cs  will  begin  to  leak  out  through  the  tube  but 
will  condense  at  the  end  farthest  from  the  side  arm  and  clog  the  tube. 
The  only  escape  route  for  the  Cs  will  be  by  diffusion  across  the  3/8- 
inch  seal  between  the  sapphire  windows  and  the  polished  plate,  which 
hopefully  will  be  quite  slow.  Thermocouple  and  heater  connections  to 
the  inner  cell  will  be  made  by  vacuum  feedthroughs  in  the  base  plate. 
Construction  of  the  new  cell  is  approaching  completion;  the  only  re- 
maining work  is  the  polishing  and  electron-beam  welding  of  the  stain- 
less-steel plate. 

We  recently  learned  of  two  different  methods  for  contacting 
sapphire  to  niobium — active  brazes  such  as  nickel-niobium  and  titanium- 
vanadium  alloys  and  calcium  aluminate  frit.  Both  the  frit  and  the 
brazes  are  extremely  resistant  to  attack  by  Cs  vapor.  Because  niobium 
is  not  attacked  by  Cs  and  its  coefficient  of  thermal  expansion  is  com- 
patible with  sapphire,  construction  of  an  all-niobium  cell  with  sapphire 
window*;  is  possible.  If  the  lifetime  of  this  cell  design  proves  to  be 
inadequate,  we  shall  pursue  the  development  of  a niobium-sapphire  cell 
that  could  be  mounted  in  the  outer  cell.  Evacuation  of  the  outer  cell 
wouli  be  required  because  niobium  is  susceptible  to  severe  oxidation  at 
elevated  temperatures . 

bong-lifetime  cells  are  particularly  important  because  a major 
goal  will  be  to  make  a careful  measurement  and  characterization  of  image 
quality.  The  image  resolution,  contrast,  coma,  and  other  aberrations 
will  be  measured  and  compared  to  theory;  methods  of  correction  where 
necessary  will  be  examined. 
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OUTSIDE  PUBLICATIONS 


1 . Papers  Presented  at  Meetings 


Information  Systems 


Cover,  T.  and  J.  van  Campenhout,  "On  the  Possible  Orderings  in  the 
Measurement  Selection  Problem,"  Third  International  Joint  Conf.  on 
Pattern  Recognition,  Coronado,  Calif.,  8-11  Nov  1976. 

Cover,  T.,  "Descriptive  Structure  of  Pattern  Recognition,"  Third 
International  Conf.  on  Pattern  Recognition,  Coronado,  Calif.,  8-11 
Nov  1976. 

Friedlander,  B.,  M.  Morf,  T.  Kailath,  and  L.  Ljung,  "Levinson-  and 
Chandrasekhar-Type  Equations  for  a General  Discrete-Time  Linear 
Estimation  Problem,"  1976  Decision  and  Control  Conf.,  Fla.,  Dec 
1976. 

Gill,  J.  and  I.  Simon,  "ink,  Dirty-Tape  Turing  Machines,  and  Quasi- 
complexity Measures,"  Third  International  Colloquium  on  Automata, 
Languages,  and  Programming,  Edinburgh,  Scotland,  Jul  1976. 

Heilman,  M.,  "Low  Cost  Data  Compression,"  International  Telemeter- 
ing Conf.,  Los  Angeles,  28-30  Sep  1976. 

Kailath,  T.,  L.  Ljung,  and  M.  Morf,  "Recursive  Input-Output  and 
State-Space  Solutions  for  Continuous -Time  Linear  Estimation  Prob- 
lems," 1976  Decision  and  Control  Conf.,  Fla.,  Dec  1976. 

Kung,  S-Y.,  T.  Kailath,  and  M.  Morf,  "A  Generalized  Resultant  Ma- 
trix for  Polynomial  Matrices,"  1976  Decision  and  Control  Conf., 
Fla.,  Dec  1976. 

Ljung,  L.  and  T.  Kailath,  "Formulas  for  Efficient  Change  of  Initial 
Conditions  in  Linear  Least  Squares  Estimation,"  1976  Decision  and 
Control  Conf.,  Fla.,  Dec  1976. 

Macovski,  A.  and  D.  Rosenfeld , "Source  Dependent  Imaging  in  Nuclear 
Medicine,"  IEEE  Symposium  on  Computer-Aided  Diagnosis  of  Medical 
Images,  Coronado,  Calif.,  11  Nov  1976. 

Morf,  M.,  L.  Ljung,  and  T.  Kailath,  "Fast  Algorithms  for  Recursive 
Identification,"  1976  Decision  and  Control  Conf.,  Fla.,  Dec  1976. 

Rosenfeld,  D.  and  A.  Macovski,  "A  New  Method  for  Three-Dimensional 
Imaging  in  Nuclear  Medicine,"  Annual  Conf.  on  Engineering  in  Medi- 
cine and  Biology,  Boston,  Mass.,  8 Nov  1976. 

57 


l 


PFECEDI1G  PAGE/BLANK-NGT  F ILM 


Digital  Systems 


Owicki,  S.  S.,  "Verifying  Monitors,"  IFIP  Working  Group  2.2,  Tamp- 
ere, Finland,  23-27  Aug  1976. 

Terman,  F.  W.,  "A  Study  of  Interleaved  Memory  Systems  by  Trace 
Driven  Simulation,"  Fourth  Symposium  on  the  Simulation  of  Computer 
Systems,  Boulder,  Colo.,  10  Aug  1976. 


Integrated  Circuits 


Barth,  P.  W.  and  J.  B.  Angell  , "Dielectrically  Isolated  Integrated 
Circuits  Made  with  Preferential  Etching,"  IEEE  International  Solid- 
State  Circuits  Conf.,  18  Feb  1977. 


Estreich,  D.  B.  and  R.  W.  Dutton,  "Modeling  Integrated  Injection 
Logic  (I-L),"  Tenth  Annual  Asilomar  Conf.  on  Circuits,  Systems,  and 
Computers,  Asilomar,  Calif.,  22  Nov  1976. 


Estreich,  D.  B.  and  R.  W.  Dutton,  "Modeling  Integrated  Injection 
Logic  (I“L),  Performance  and  Operational  Limits,"  International 
Solid-State  Circuits  Conf.,  Philadelphia,  16  Feb  1977. 


Solid  State 


Garner,  C.  M.,  Y.  D.  Shen , G.  L.  Pearson,  and  W.  E.  Spicer,  "Auger 
Profiling  of  Interface  Widths  of  AlxGa^_xAs  Hetero junctions ,"  meet- 
ing of  the  Am.  Phys.  Soc.,  Stanford  University,  Stanford,  Calif., 
20-22  Dec  1976.  Abstract  published  in  Bull.  Am.  Phys.  Soc.,  21, 
1976,  p.  1290. 

Garner,  C.  M.,  Y.  D.  Shen,  G.  L.  Pearson,  and  W.  E.  Spicer,  "Auger 
Profiling  of  Interface  Widths  of  AlxGa^_xAs  Hetero junctions  Grown 
by  LPE,"  meeting  of  the  Am.  Phys.  Soc.,  Stanford  University,  Stan- 
ford, Calif.,  20-22  Dec  1976. 

Houng,  Y.  M.,  B.  L.  Mattes,  and  G.  L.  Pearson,  "Deep  Trapping  Ef- 
fects at  the  GaAs-GaAs:Cr  Interface  in  GaAs  FET  Structures,"  meet- 
ing of  the  Am.  Phys.  Soc.,  Stanford  University,  Stanford,  Calif., 
20-22  Dec  1976. 
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This  appendix  is  a compilation  of  abstracts  of  reports  issued  by 
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A STUDY  OF  INTERLEAVED 
MEMORY  SYSTEMS  BY 
TRACE  DRIVEN  SIMULATION 
by  Fred  W,  Terman 


JSEP  Contract 
N00014-75-C-0601 
Tech  . Note  No  . 94 
September  1976 


ABSTRACT 


A model  of  interleaved  memory  systems  for  IBM  360/370  architecture 
has  been  investigated  by  means  of  a trace  driven  simulation.  The  model 
used  is  an  extension  of  one  due  to  G.  J.  Burnett  and  E.  G.  Coffman,  Jr. 
The  trace  data  to  drive  the  simulation  was  obtained  from  traces  of  typ- 
ical IBM  360/370  programs  and  of  the  0S/VS2  operating  system.  The  pre- 
dictions of  the  Burnett -Coffman  model  are  found  to  fit  well  with  the 
simulation  results  for  the  fetching  of  instructions.  For  the  fetching 
of  operands,  however,  the  simulation  results  show  only  about  half  the 
increase  in  the  memory  bandwidth  that  the  Burnett-Cof fman  analysis  pre- 
dicts. This  indicates  that  data  references  on  the  IBM  360/370  are  not 
random  in  the  way  assumed  by  Burnett  and  Coffman. 
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MAGNETOS PHKRIC  DUCT 

PROPAGATION 

by  Paul  A.  Bernhardt 


N00014-76-C-0689 
TR  No.  3477-1 
October  1976 


ABSTRACT 

New  propagation  models  have  been  developed  for  the  study  of  the 
ducting  of  whistlers.  Two  classes  of  models  are  used  in  this  study.  One 
class  of  models  is  based  on  the  electromagnetic  laws  governing  the  wave 
propagation.  The  other  class  describes  the  protonospheric  and  iono- 
spheric plasma  concentrations  which  affect  the  trajectories  of  VLF  waves. 
Since  the  electromagnetic  laws  governing  VLF  propagation  have  been  more 
fully  explored  than  the  physical  laws  governing  the  environment  in  which 
VLF  waves  propagate,  we  have  concentrated  on  developing  new,  realistic 
models  of  the  plasmasphere  and  have  used  standard  , raytracing  techniques 
to  investigate  magnetos pheric  duct  propagation. 

The  first  part  of  this  report  is  a description  of  a protonospher ic- 
ionospheric  model  used  to  simulate  the  diurnal  and  seasonal  behavior  of 
the  magnetosphere.  Special  attention  is  given  to  the  protonospheric- 
ionospheric  interface.  This  is  the  region  where  enhancement  ducts, 
which  trap  VLF  waves  in  the  magnetosphere,  merge  into  the  relatively 
uniform  ionosphere.  The  altitude  where  this  merging  occurs  defines  the 
duct  end-point  height.  Temporal  variations  in  the  height  of  the  end- 
points of  a magnetnspheric  duct  affect  the  coupling  between  the  earth's 
surface  and  the  duct  and  the  echoing  properties  of  the  duct.  The  rela- 
tionship between  0+ - H+  ion  transition  height  and  the  end-points  of  mag- 
netospheric  ducts  is  described.  The  decay  of  a whistler  duct  due  to 
thermal  plasma  flow  is  also  discussed. 

The  second  part  of  this  report  is  a study  of  wave  trapping  proper- 
ties of  magnetospheric  ducts.  Rays  are  injected  at  the  geomagnetic 
equator  with  a range  in  wave  normal  angles  and  positions  inside  a duct. 
The  altitude  where  the  rays  leave  the  duct  is  a measure  of  the  wave 
trapping  efficiency. 
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FINAL  REPORT  ON 
COMPUTER-AIDED  SEMICON- 
DUCTOR PROCESS  MODELING 
by  J.  D.  Meindl , Director 


ABSTRACT 


University  of  Florida 
Contract  DAAB07-75-C-1344 
from  the  ARPA 
1 Jul  1975  - 30  Sep  1976 


In  military,  industrial,  commercial,  and  consumer  applications , 
frequently  there  is  a great  necessity  for  "customizing”  the  design  of 
an  integrated  circuit  to  fulfill  the  critical  needs  of  a specific  sys- 
tem or  class  of  systems.  A major  barrier  which  prevents  the  economic 
production  of  small  quantities  of  high  performance  custom  integrated 
circuits  is  the  cost  of  design.  That  is,  the  initial  cost  of  designing 
optimum  fabrication  processes,  device  structures,  and  circuit  configu- 
rations is  prohibitively  large  because  of  the  amount  of  empirical  human 
effort  which  must  be  invested.  The  root  cause  of  this  problem  is  a 
glaring  lack  of  adequate  process,  device,  and  circuit  models  and  accom- 
panying computer  aided  design  techniques  to  ease  the  burden  of  custom 
design.  Perhaps  the  most  serious  deficiency  among  these  is  the  unavail- 
ability >f  suitable  models  for  predicting  accurately  the  physical  char- 
acteristics of  a monolithic  structure  on  the  basis  of  the  control  param- 
eters for  the  correspond ing  fabrication  processes.  The  salient  objec- 
tive of  this  program  is  the  development  of  new  basic  models  for  inte- 
grated circuit  processes  which  will  permit  accurate  prediction  of  the 
characteristics  of  a monolithic  structure  on  the  basis  of  its  proposed 
process  parameters.  These  models  will  serve  as  the  basis  for  economic 
computer  aided  design  of  optimum  fabrication  processes  for  custom  inte- 
grated circuits. 

The  four  key  generic  integrated  circuit  fabrication  processes  which 
are  being  investigated  are  (1)  ion  implantation,  (2)  thermal  oxidation 
and  chemical  vapor  deposition,  (3)  epitaxy,  and  (4)  thermal  diffusion. 
Tills  report  describes  the  progress  which  has  been  made  in  the  first  15 
months  of  this  program  in  all  four  areas. 

In  ion  implantation,  studies  have  been  done  in  two  areas.  Effects 
of  implantation  damage  on  impurity  profiles  in  annealed  Si  have  been 
investigated.  A three-stream  diffusion  model  for  boron  in  silicon  in- 
volving boron,  vacancies  and  boron  vacancy  complexes  has  been  proposed. 
The  model  is  capable  of  predicting  ordinary  diffusion,  proton  enhanced 
diffusion  and  the  annealing  behavior  of  room-temperature  implanted  boron 
when  appropriate  restrictions  on  dose  and  annealing  temperature  are 
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obeyed.  In  the  second  study  a straightforward  approach  to  calculations 
of  range  profiles  of  dopant  ions  in  multilayer  media  (e.g.  arsenic  im- 
plementation in  silicon  through  a thin  layer  of  silicon  dioxide)  has 
been  developed  based  on  numerical  integration  of  the  linearized  Boltz- 
mann transport  equation.  Using  this  approach,  we  have  also  calculated 
the  range  profiles  of  the  recoiling  ions  which  are  knocked  on  from  one 
layer  of  the  substrate  into  the  other  layer  by  the  projectiles  (e.g. 
oxygen  recoiling  from  SiO^  into  Si,  by  ions). 

In  the  area  of  thermal  oxidation  and  chemical  vapor  deposition  the 
principal  goal  of  the  first  year  of  the  program  has  been  to  achieve  ac- 
curate analytic  prediction  of  oxide  thickness  for  an  arbitrary  sequence 
of  oxidations.  At  the  present  time,  such  prediction  is  possible  only 
for  lightly  doped,  (111)  oriented  silicon,  with  no  chlorine  species 
present  during  the  oxidation  process. 

Toward  this,  kinetic  oxidation  data  have  been  gathered  under  the 
following  conditions: 

(a)  Data  for  (111)  and  (100)  silicon  over  the  temperature 
range  700 °C  to  1200°C. 

(b)  Data  for  dry  O2/HCL  oxidation  ambients  containing  0-10*^ 

HCL  for  temperatures  between  900  °C  and  1100 °C . 

(c)  Data  for  heavily  phosphorous  doped  (111)  substrates  for 
surface  concentrations  up  to  solid  solubility  and  tem- 
peratures between  900 °C  and  1100 °C. 

The  availability  of  the  above  data  now  makes  it  possible  to  predict 
analytically  the  oxide  thickness  under  a very  wide  range  of  ambient  con- 
ditions and  substrate  orientation  and  doping  levels. 

In  epitaxy,  preliminary  effort  has  been  directed  toward  understand- 
ing the  kinetics  of  growth  and  the  mechanism  of  dopant  incorporation  in 
the  epitaxial  layer. 

The  dopant  system  of  a horizontal  silicon  epitaxial  reactor  has 
been  characterized  by  changing  the  dopant  gas  flow  during  the  continuous 
deposition  of  epitaxial  layers  from  silane.  A system  "transfer  function" 
has  been  found  relating  the  dopant  profile  in  the  epitaxial  layer  (system 
output)  to  the  time-varying  dopant  gas  flow  (system  input).  The  transfer 
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function  allows  the  calculation  of  the  dopant  profile  in  the  epitaxial 
film  for  any  time  varying  dopant  gas  flow.  The  calculation  of  the  dop- 
ant gas  flow,  as  a function  of  the  time  required  to  achieve  a desired 
dopant  profile,  is  thus  possible. 

In  thermal  diffusion,  a new  mathematical  model  is  being  developed 
for  the  diffusion  of  impurities  into  silicon  for  both  a constant  source 
and  drivein  diffusion  process.  Tbe  model  considers  the  influence  of 
the  internal  electric  field  on  the  motion  of  impurity  ions  at  elevated 
temperatures.  The  electric  field  is  due  to  charge  density  produced  by 
the  ionization  of  impurities  present  in  the  material.  Appendix  B re- 
views the  results  of  this  study. 
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ERROR  CORRECTION  BY 
A I.TERNATE-DATA  RETRY 
by  John  J.  Shedletsky 


IBM  Grad.  Fellowship 
NSF  Grant  GJ-40286 
TR  No.  113 
May  1976 


ABSTRACT 

A new  technique  for  low-cost  error  correction  in  computers  is  the 
alternate-data  retry  (ADR).  An  ADR  is  initiated  by  the  detection  of  an 
error  in  the  initial  execution  of  an  operation.  The  ADR  is  a re-execu- 
tion  of  the  operation,  but  with  an  alternate  representation  of  the  ini- 
tial data.  The  choice  of  the  alternate  representation  and  the  design 
of  the  processing  circuits  combine  to  insure  that  even  an  error  due  to 
a permanent  fault  is  not  repeated  during  retry.  Error-correction  is 
provided  at  a hardware  cost  comparable  to  that  of  a conventional  retry 
capability. 

Sufficient  conditions  are  given  for  the  design  of  circuits  with  an 
ADR  capability.  The  application  of  an  ADR  capability  to  memory  and  to 
the  data  paths  of  a processor  is  illustrated. 
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CASCADE  STRUCTURE  NSF  Grant  GJ-27527 
IN  TOTALLY  SELF-  TR  No.  108 
CHECKING  NETWORKS  April  1976 
by  Stephen  Kolupaev 


ABSTRACT 


In  the  well-known  totally  self-checking  (TSC)  network,  a failure 
must  not  change  one  output  codework  into  another.  Called  the  fault-se- 
cure property,  this  permits  a receiver  of  the  net's  output  to  assume 
that  any  codework  it  receives  is  correct.  Further,  the  self-testing 
property  requires  that  each  possible  failure  in  the  net  must  produce  at 
least  one  non-code  output.  Thus  a receiver  can  monitor  the  health  of 
the  network  by  watching  for  non-code  outputs. 

In  this  paper  we  propose  modifications  of  these  two  properties.  The 
self-testing  property  is  made  more  stringent.  Each  possible  failure  in 
the  net  is  required  to  produce  an  output  which  is  in  a distinguished  sub- 
set  of  the  non-code  outputs.  The  fault-secure  requirement  is  modified 
to  permit  a fault  to  interchange  certain  output  codewords.  In  particu- 
lar, all  outputs  not  in  the  distinguished  subset  are  partitioned  into 
equivalent  classes,  and  a fault  is  permitted  to  change  the  output  from 
one  codeword  to  another  codeword  in  the  same  class.  However,  a fault  is 
not  permitted  to  change  the  output  from  a codeword  to  any  member  of  a 
different  equivalence  class  (one  not  containing  the  correct  output). 

These  modified  properties  define  a generalization  of  the  TSC  net- 
work. A network  which  meets  the  modified  properties  is  called  a gener- 
alized self-checking  (GSC)  network.  Self-checking  and  self-testing 
(Morphic)  networks  and  TSC  networks  are  special  cases  of  the  GSC  network. 

Examining  TSC  networks,  we  find  a further  connection  with  the  GSC 
network.  It  has  been  known  for  some  time  that  not  every  subnetwork  of  a 
TSC  networ’/  need  be  TSC.  We  show  that  every  subnetwork  of  a TSC  network 
is  GSC,  and  every  TSC  network  is  a cascade  of  GSC  networks.  This  estab- 
lishes the  GSC  network  as  the  basic  building  block  from  which  every  TSC 
network  is  constructed . 

We  explore  a brute-force  method  for  constructing  a desired  TSC  net- 
work by  cascading  GSC  subnetworks.  The  method  resorts  to  enumeration  at 
many  points  of  decision  and  thus  is  not  a practical  design  tool.  However, 
it  does  yield  a very  nice  alternate  realization  of  the  Morphic  OR,  and 
suggests  specializations  which  merit  further  study. 
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SYMMETRY,  AUTOMORPHISM, 
AND  TEST 

by  Andre  Verdillon 


Inst,  de  Recherche  en 
Inf.  et  en  Auto.  Grant 
NSF  Grant  MCS  76-05327 
Tech.  Note  No.  87 
June  1976 

ABSTRACT 

The  study  of  regular  structures  such  as  tree,  Iterative,  monotonic, 
functions  has  shown  that  by  taking  into  account  their  properties  the 
problem  of  finding  a test  can  be  reduced  considerably.  This  paper  shows 
how  testing  can  be  simplified  by  taking  account  symmetries  that  lead  to 
network  automorphisms.  For  an  acyclic  network  of  gates,  methods  are 
presented  that  give  either  an  optimal  or  a near  optimal  test.  It  is 
shown  how  the  concept  of  automorphism  can  be  extended  to  networks  of 
modules  and  how  this  allows  using  a multi-level  description  of  a system. 
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ERROR  MANAGEMENT 
IN  DIGITAL  COMPUTER 
INPUT/OUTPUT  SYSTEMS 
by  Alan  Michael  Usas 


NSF  Fellow 
NSF  Grants 
GJ— 10286  and  GK-43322 
TR  No.  122 
May  1976 


ABSTRACT 


Historically  peripheral  devices  have  been  more  important  in  deter- 
mining the  market  success  of  a computer  system  than  the  central  proces- 
sor. In  spite  of  their  key  role  the  input/output  (I/O)  system  components 
have  been  encumbered  with  a poor  reputation  for  reliability  and  maintain- 
ability. To  date  there  has  been  no  comprehensive  treatment  of  this  prob- 
lem. Presented  here  are  techniques  of  error  management  (those  steps 
taken  to  limit  the  effect  of  errors  and  to  expedite  recovery)  appropriate 
for  the  buses,  interfaces,  controllers,  and  peripheral  devices  of  the  I/O 
system . 

Many  of  the  possible  erroneous  operations  in  an  I/O  system  can  be 
classified  as  either  intolerable  or  tolerable  on  the  basis  of  whether 
their  effect  is  to  physically  damage  a part  of  the  system,  such  as  a 
head  "crash"  on  a disk  drive,  or  cause  some  other  irreversible  action, 
or  whether  the  error  simply  renders  a device  inoperable  without  destroy- 
ing any  information.  Fail-safe  circuits  are  designed  to  minimize  the 
probability  of  just  one  of  the  two  possible  erroneous  outputs  from  a 
single  output  circuit  (i.e.,  erroneous  1 and  0),  and  thus  represent  a 
useful  compromise  solution  over  circuits  that  would  reduce  the  likeli- 
hood of  both  types  of  errors.  The  design  of  fail-safe  combinatorial 
circuits  by  use  of  asymmetrically  failing  components  and  circuit  dupli- 
cation 1 explored,  and  a new  synthesis  algorithm  is  presented.  It  is 
based  on  the  probability  of  an  unsafe  erroneous  output,  which  is  a func- 
tion of  the  input  vectors  and  possible  faults,  rather  than  on  just  the 
probability  of  component  failure.  Fail-safe  networks  constructed  using 
this  algorithm  will  generally  require  less  than  complete  duplication  and 
will  have  a lower  probability  of  intolerable  erroneous  outputs  than  if 
fault  probabilities  alone  were  considered.  As  sequential  circuits  are 
found  throughout  the  interfaces  and  controllers  of  an  I/O  system,  meth- 
ods of  constructing  fail-safe  versions  that  will  not  produce  intolerable 
errors  even  after  an  erroneous  state  has  been  entered  are  explored . A 
straightforward  design  requiring  at  most  two  redundant  state  variables 
that  explicitly  traps  the  circuit  in  an  erroneous  state  identified  by 
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means  of  an  error-detecting  code  is  described  and  shown  to  be  suitable 
for  large-scale  integration. 

Besides  the  direct  prevention  of  intolerable  errors,  the  problem 
of  detecting  the  occurrence  of  safe  but  nevertheless  erroneous  condi- 
tions is  studied.  To  prevent  checker  or  error  detector  failures  from 
masking  the  indication  of  an  erroneous  condition,  these  circuits  should 
fail  to  the  error-indicating  state.  Designs  for  a variety  of  fail-safe 
checkers  for  velocity,  position,  timing,  and  other  important  peripheral 
device  parameters  that  utilize  the  failure  asymmetry  of  electro-optical 
components  or  duplication  are  described. 

Specific  modifications  of  magnetic  disks,  tapes,  card  readers  and 
punches,  line  printers,  and  terminals  using  fail-safe  circuits  to  re- 
duce the  probability  of  an  intolerable  error  and  checkers  to  enable 
timely  error  detection  are  detailed.  A similar  treatment  is  given  to 
device  controllers  and  interfaces  with  the  further  purpose  of  limiting 
the  propagation  of  errors  and  that  utilizes  the  controller  to  verify 
proper  peripheral  device  operation. 

A flexible  maintenance  bus  is  suggested  as  a solution  to  the  prob- 
lem of  simplifying  the  human  involvement  with  I/O  system  fault  diag- 
nosis. This  bus,  which  is  physical! v attached  to  a suspected  faulty 
piece  of  equipment,  acts  like  a programmable  waveform  generator  and 
response  recorder.  It  gives  diagnostic  programs  access  to  a multiplic- 
ity of  low  level  points  in  the  I/O  system  not  directly  reachable  by 
means  of  the  standard  I/O  bus . 
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ABSTRACT 


The  pulse-echo  Imaging  modality  employed  in  medical  acoustic  imag- 
ing systems  permits  direct  three-dimensional  visualization  of  a reflect- 
ing object  embedded  in  a three-dimensional  acoustic  medium  (such  as  an 
organ  in  the  interior  of  the  human  body) . High  resolution  in  range  is 
achieved  by  range-gating  ultrasonic  pulses  and,  for  the  more  advanced 
systems,  high  resolution  in  azimuth  is  achieved  by  focusing  the  sound 
with  the  aid  of  an  array  of  acoustic  transducers  coupled  to  a system  of 
controlled  analog  delay  elements;  the  dynamic  variation  of  these  delay 
elements  permits  systematic  exploration  with  the  focused  sound  over  the 
region  of  interest  in  the  object  space.  In  such  a fashion,  acoustic  re- 
flectivity is  mapped  into  a brightness  pattern  (the  image)  on  the  face 
of  a cathode  ray  tube. 

Existing  and  proposed  electronically-focused  schemes  of  this  kind 
are,  with  few  exceptions,  based  on  delay/sum  processing  of  the  detected 
acoustic  signals;  i .e . , the  outputs  of  the  delay  elements  are  coherently 
added  together  to  yield  one  point  in  the  Image  at  a time.  Until  recently, 
most  studies  of  focused,  pulse-echo  acoustic  imaging  systems  have  assumed 
that  a steady-state  (monochromatic)  analysis  was  adequate  to  derive  the 
essential  characteristics  of  the  system's  performance,  such  as  axial  and 
lateral  resolution  and  grating  lobe  properties.  When  very  wideband 
pulses  are  employed  (to  achieve  high  depth  resolution),  however  a more 
general  wideband,  temporal -spatial  analysis  of  the  system  point  response 
is  required.  For  example,  as  the  acoustic  pulses  become  increasingly 
wideband,  it  is  shown  in  this  research  that  the  lateral  response  of  the 
(two-dimensional)  temporal-spatial  point  spread  function  provided  by  sim- 
ple delay/sum  processing  becomes  increasingly  poor.  Here,  a method  is 
developed  for  restoring  the  degraded  lateral  response  that  takes  place 
under  very  wideband  conditions;  this  scheme  is  realized  by  a one-dimen- 
sional filtering  operation  performed  on  the  output  of  the  delay/sum  stage 
of  the  basic  imaging  system.  This  leads  to  an  improved  overall  axial- 
lateral  (temporal -spatial)  point  spread  function  for  the  electronically- 
focused  pulse-echo  imaging  system. 
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To  develop  the  background  for  the  study  of  wideband  acoustic  imaging 
systems,  the  theory  of  monochromatic  scalar  diffraction  is  extended  to 
wideband  d isturbances . A variety  of  applications  of  the  wideband  theory 
is  provided  for  illustration.  Among  these  is  a description  of  the  wide- 
band focusing  of  a plane  wave  impulse  by  a circular  lens. 

Following  this,  the  subject  of  multiple-element  focused  array  sys- 
tems is  developed.  First,  the  monochromatic  response  of  the  focused 
system  is  derived;  then,  the  theory  of  array  focusing  is  extended  to 
wideband  operation.  In  particular,  the  wideband  point  response  is  ob- 
tained, and  its  dependence  on  the  form  of  the  temporal  pulse  is  empha- 
sized. Also,  wideband  grating-lobe  behavior  is  investigated,  and  a con- 

th 

venient  relative  measure  of  comparison  between  the  strength  of  the  n 
grating  lobe  and  the  main  lobe  is  derived. 

Next,  the  general  sub.ject  of  wideband  image  reconstruction  is  de- 
veloped, where  delay/sum  processing  is  evaluated  and  the  axial-lateral 
resolution  optimizing  scheme  is  derived.  Here,  two  possible  reconstruc- 
tion strategies  are  contrasted:  coherent  summation  and  deconvolution. 

It  is  shown  that  deconvolution  provides  the  optimum  wideband  point  re- 
sponse in  the  sense  that  an  optimum  compromise  between  axial  and  lateral 
resolution  is  automatically  obtained,  in  contrast  to  coherent  summation. 
The  subject  of  wideband  image  reconstruction  is  concluded  with  the  in- 
troduction of  an  alternative  formulation  of  the  problem  of  image  recon- 
struction by  deconvolution:  an  analytical  inversion  formula  is  derived 

that  allows  a function  of  reflectivity  (the  object)  to  be  recovered  from 
measurements  of  its  line  integrals  defined  over  circular  paths. 

Finally,  the  subject  of  transmit/receive  aperture  synthesis  is  ex- 
amined. A simple  rule  is  derived  which  shows  how  a combination  of  two 
"incomplete"  or  sparse  transmitting  and  receiving  structures  can  be 
exploited  to  provide  full  aperture  resolution.  A variety  of  original 
examples  is  provided  to  illustrate  this  principle. 
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by  M.  Heilman,  R.  Merkle, 

R.  Schroeppel , L.  Washington, 

W.  Dlffie,  S.  Pohlig , 

P.  Schweitzer 


ABSTRACT 

The  National  Bureau  of  Standards  has  proposed  a data  encryption 
standard  (DES)  for  nonmilitary  governmental  and  commercial  use.  Such  a 
cryptosystem  can  be  used  to  provide  both  privacy  and  authentication  to 
messages  [1,2]  . In  many  applications,  authentication  (a  form  of  digital 
signature)  is  of  primary  interest. 

For  reasons  of  national  security,  NBS,  IBM  (which  designed  DES), 
and  the  National  Security  Agency  (which  participated  in  its  evaluation), 
will  not  divulge  the  results  of  their  tests  concerning  the  strength  of 
the  system.  Because  they  will  also  not  explain  the  rationale  behind 
certain  critical  elements  or  adequately  quantify  the  level  of  security, 
we  undertook  an  evaluation  of  DES  in  an  attempt  to  determine  its 
strengths  and  weaknesses.  This  evaluation  lasted  one  month  and  involved 
approximately  ten  person  weeks  of  effort.  Our  results  are  therefore 
only  preliminary,  and  it  must  be  assumed  that  an  adversary  with  greater 
resources  would  make  greater  progress . 

It  is  thus  somewhat  disturbing  that  we  achieved  a 50  percent  sav- 
ings in  computation  over  an  exhaustive  search.  We  also  found  suspicious 
structure  in  the  critical  elements  of  the  system  (the  S-boxes).  This 
structure  could  be  the  result  of  an  accidental  weakness,  a deliberately 
set  trap  door  (i.e.,  hidden  structure  which  allows  those  who  know  of  it 
to  break  the  system),  or  no  weakness  at  all.  An  explanation  and  further 
study  are  needed  before  trust  can  be  placed  in  DES.  This  need  is  en- 
hanced because  NSA  does  not  want  a genuinely  strong  system  to  frustrate 
its  cryptanalytic  intelligence  operations.  As  a result,  DES  is  mildly 
suspect  a priori. 
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ABSTRACT 

Reliability  is  an  important  concern  of  any  computer  system.  No 
matter  how  carefully  designed  and  constructed,  computer  systems  fail. 

The  rapid  and  systematic  restoration  of  service  after  an  error  or  mal- 
function is  always  a major  design  and  operational  goal.  In  order  to 
overcome  the  effects  of  a failure,  recovery  must  be  performed  to  go 
from  the  failed  state  to  an  operational  state.  This  thesis  describes 
a recovery  method  which  guarantees  that  a computer  system,  its  associ- 
ated data  bases  and  communication  transactions  will  be  restored  to  an 
operational  and  consistent  state  within  a given  time  and  cost  bound 
after  the  occurrence  of  a system  failure. 

This  thesis  considers  the  optimization  of  a specific  software 
strategy  - the  rollback  and  recovery  strategy,  within  the  framework  of 
a graph  model  of  program  flow  which  encompasses  communication  interfaces 
and  data  base  transactions.  Algorithms  are  developed  which  optimize 
the  placement  of  dynamic  recovery  checkpoints . Presented  is  a method 
for  statically  pre-computing  a set  of  optimal  decision  parameters  for 
the  associated  program  model,  and  a run-time  technique  for  dynamically 
determining  the  optimal  placement  of  program  recovery  checkpoints  . 
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EQUIDENSITOMETRY  BY 
COHERENT  OPTICAL 
FILTERING 

by  H.  K.  Liu,  J.  W. 
Goodman,  J.  L-H . Chan 


ABSTRACT 


A new  and  simple  method  has  been  developed  for  fabrication  of  mul- 
tilevel halftone  screens  that  have  proved  useful  for  generation  of  equi- 
density  contours  by  means  of  coherent  optical  filtering.  Up  to  ten  con- 
tours of  constant  brightness  have  been  generated  on  a picture  by  coherent 
optical  filtering  of  a single  halftone  copy  of  that  picture.  The  possible 
application  of  this  method  to  optical  analog-to-digital  conversion  using 
a single  halftone  photograph  is  discussed. 


iii 


SEL-76-044 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Oata  Enteredl 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


2.  GOVT  ACCESSION  NO  I 3 RECIPIENT  S CATALOG  NUMBER 


REPORT  DOCUMENTATION  PAGE 


1 REPORT  NUMBER 


5 TYPE  OF  REPORT  N PERIOD  COVERED 

Semi-Annual  Status  Report 
1 Jul  through  31  Dec  1976 


6 PERFORMING  ORG  REPORT  NUMBER 

7 author..,  SEL-76-045 ^ 

. , a.  CONTRACT  OR  GRANT  NUMBER!.) 

Information  Systems,  Digital  Systems,  Integrated  . £- 

Circuits,  Solid-State  Electronics,  Radioscience,  ^ ~^^00014-75-C-0601 

Hansen  Laboratories,  and  Inst,  for  Plasma  Research  — J — - 


9 PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Stanford  Electronics  laboratories 
Stanford  University  ^ - 

Stanford,  California  94305 


11  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

U.S.  Navy 

Office  of  Naval  Research 


14  MONITORING  AGENCY  NAME  8,  ADDRESS  (if  diff.  from  Controlling  0„ic«) 


10  PROGRAM  ELEMENT,  PROJECT  TASK 
AREA  & WQR.K_UKJI  NUMBERS 

1 1 i Dec-  y c 


12  REPORT  DF)T£  -19  NO  OF  PAGES 

None  78 


15  SECURITY  CLASS  (of  thj.  'ejMM4f ■*.. 

Unclassified  ^ ^ fO^LJ  I 


15a  DECLASSIF  ICATION /DOWNGRADING 
SCHEDULE 


16  DISTRIBUTION  STATEMENT  lof  thi.  report) 

This  document  has  been  approved  for  public  release  and  sale;  its  distribution 
is  unlimited.  Reproduction  in  whole  or  in  part  is  permitted  for  any  purpose 
of  the  United  States  Government . 


17  DISTRIBUTION  STATEMENT  (of  the  abstract  antarad  in  Block  20.  if  diffarant  from  report) 


19  KEY  WOPDS  (Contmua  on  ravarta  tide  if  necettary  and  idantify  by  block  number) 


PLASMA  PHYSICS,  QUANTUM  ELECTRONICS,  RADIOSCIENCE,  SOLID-STATE  ELECTRONICS, 
SYSTEMS  THEORY,  DIGITAL  SYSTEMS,  INFORMATION  SYSTEMS 


20  ABSTRACT  (Contmua  on  ravarta  tide  if  nacataary  and  identify  by  block  number) 

The  progress  of  research  under  each  of  the  locally  established  project  numbers 
for  this  contract  is  summarized . 


DD,'0“„1473 

EDITION  OF  1 NOV  «S  IS  OBSOLETE 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  IWhpn  Data  Entenwll 


JSEP  REPORTS  DISTRIBUTION  LIST 


No.  of 
Copies 


Department  of  Defense 

Defense  Documentation  Center 
Attn:  DDC-TGA  (Mrs.  V.Caponio)  12 

Cameron  Station 
Alexandria,  Virginia  22314 

Asst.  Dir. , Electronics  1 

and  Computer  Sciences 
Office  of  Director  of  Defense 
Research  and  Engineering 
The  Pentagon 
Washington,  D. C.  20315 

Office  of  Director  of  Defense 
Research  and  Engineering 
Information  Office  Lib.  Branch 
The  Pentagon 
Washington,  D.C.  20301 

ODDR&E  Advisory  Group  on 
Electron  Devices 
201  Varick  Street 
New  York,  New  York  10014 

Chief,  R&D  Division  (340) 

Defense  Dommunications  Agency 
Washington,  D.C.  20301 

Director,  Nat.  Security  Agency 
Fort  George  G.  Meade 
Maryland  20755 
Attn:  Dr.  T. J. Beahn 

Institute  for  Defense  Analysis 
Science  and  Technology  Division 
400  Army-Navy  Drive 
Arlington,  Virginia  22202 

Dr.  Stickley 

Defense  Advanced  Research 
Projects  Agency 
Attn:  Technical  Library 

1400  Wilson  Boulevard 
Arlington,  Virginia  22209 


Dr.  R.  Reynolds 
Defense  Advanced  Research 
Projects  Agency 
Attn:  Technical  Library 

1400  Wilson  Boulevard 
Arlington,  Virginia  22209 

Department  of  the  Air  Force 

AF/RDPS 
The  Pentagon 
Washington,  D.C.  20330 

1 AFSC  (LJ/  Mr.  Irving  R.  Mirman) 

Andrews  Air  Force  Base 
Washington,  D.C.  20334 

Directorate  of  Electronics 
and  Weapons 
1 HQ  AFSC/DLC 

Andrews  AFB,  Maryland  20334 

Directorate  of  Science 
HQ  AFSC/DLS 

1 Andrews  Air  Force  Base 

Washington,  D.C.  20334 

LTC  J. W.  Gregory 
1 AF  Member,  TAC 

Air  Force  Office  of 
Scientific  Research 
Bolling  Air  Force  Base 
Washington,  D.C.  20332 

1 

Mr.  Carl  Sletteu 
RADC/ETE 

Hanscom  AFB,  Maryland  01731 

1 Dr.  Richard  Picard 

RADC/ETSL 

Hanscom  AFB,  Maryland  01731 

Mr.  Robert  Barrett 
RADC/ETS 

Hanscom  AFB,  Maryland  01731 


-1- 

JSEP  3/77 


No.  of 
Copies 

1 


1 


1 


1 


1 


5 


1 


1 


1 


No.  of 
Copies 


No.  of 
Cop i es 


Dr.  John  N.  Howard  1 

AFGL/CA 

Honscom  AFB,  Maryland  01731 

Dr.  Richard  B.  Mack  1 

RADC/ETER 

Hanscom  AFB,  Maryland  01731 

Documents  Library  (TILD)  1 

Rome  Air  Development  Center 
Griffiss  AFB,  New  York  13441 

Mr.  H. E. Webb,  Jr.  (ISCP)  1 

Rome  Air  Development  center 
Griffiss  AFB,  New  York  13441 

Mr.  Murray  Kesselman  (ISCA)  1 

Rome  Air  Development  Center 
Griffiss  AFB,  New  York  13441 

Mr.  W.  Edwards  1 

AFAL/TE 

Wright-Patterson  AFB 
Ohio  45433 


Mr.  R,  D.  Larson  1 

AFAL/DHR 

Wright-Patterson  AFB 
Ohio  45433 

Howard  H. Steenbergen  1 

AFAL/DHE 

Wright-Patterson  AFB 
Ohio  45433 

Chief  Scientist  1 

AFAL/CA 

Wright-Patterson  AFB 
Ohio  45433 

HQ  ESD  (DRI/Stop22)  1 

Hanscom  AFB,  Maryland  01731 

Professor  R.E.  Fontana  1 

Head,  Dept,  of  Electrical  Engr. 
AFIT/ENE 

Wright-Patterson  AFB 
Ohio  45433 


Mr.  John  Mottsmith  (MCIT)  1 

HQ  ESD  (AFSC ) 

Hanscom  AFB,  Maryland  01731 

L'PC  Richard  J.  Gowen  1 

Professor 

Dept,  of  Electrical  Engineering 
USAF  Academy,  Colorado  80840 


AUL/LSE-9663  1 

Maxwell  AFB,  Alabama  36112 

AFETR  Technical  Library  1 

P.0.  Box  4608,  MU  5650 
Patrick  AFB,  Florida  32542 

ADTC  (DLOSL)  1 

Eglin  AF’B,  Florida  32542 

HQ  AMD  ( KDR/Col . Godden ) 1 

Brooks  AFB,  Texas  78235 

USAF  European  Office  of  1 

Aerospace  Research 


Technical  Information  Office 
Box  14,  FPO,  New  York  09510 

Dr.  Carl  E.  Baum  1 

AFWL  (ES) 

Kirtland  AFB,  New  Mexico  87117 

ASAFSAM/RAL  1 

Brooks  AFB,  Texas 

Department  of  the  Army 

HQ  DA  (DAMA0ARZ-A)  1 

Washington,  D.C.  20310 

Commander  1 

U. S.  Army  Security  Agency 
Attn:  IARD-T 

Arlington  Hall  Station 
Arlington,  Virginia  22212 


-2- 


JSEP  3/77 


No.  of 
Copies 

Commander  U.  S.  Army  Materiel  1 

L)ev.  & Readiness  Command 
Attn:  Tech.  Library  Rm  7S  35 

5001  Eisenhower  Ave. 

Alexandria,  Virginia  22333 

Commander  Research  Laboratory  1 

ATTN.  DRXRD-BAD 
U.S.  Army  Ballistics 
Aberdeen  Proving  Ground 
Aberdeen,  Maryland  21005 

Commander  1 

Picatinny  Arsenal 
Dover,  New  Jersey  07081 

ATTN:  SMUPA-TS-T-S 

ATTN:  Dr.  Herman  Robl  1 

U.S.  Army  Research  Office 
P. 0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 

ATTN:  MR.  Richard  0.  Ulsh  1 

U.S.  Army  Research  Office 
P.0.  Box  12211 
Research  Triangle  Park 
North  Carolina  27709 

Mr.  George  C. White,  Jr.  1 

Deputy  Director 
Pitman-Dunn  laboratory 
Frankford  Arsenal 
Philadelphia,  Penna.  19137 

Commander  1 

Attn:  Chief,  Document  Section 
U.S.  Army  Missile  Command 
Redstone  Arsenal,  Alabama  35809 

Commander  1 

U.S.  Army  Missile  Command 
Attn:  DRSMI-RR 

Redstone  Arsenal,  Alabama  35809 
C • immander 

' ■ i *•  f . Materials  Sciences 
• » l ■ «n.  Bldg.  292 

Materials  and  Mechanics 
«•  ■ ut'  ' Center 

« ••  Massachusetts  02172 


No.  of 
Copies 

Commamder 

Harry  Diamond  Laboratories  1 

ATTN:  Mr.  John  E.  Rosenberg 

2800  Posder  Mill  Road 
Adelphi , Maryland  20783 

Commandant 

U.S.  Army  Air  Defense  School  1 

Attn:  ATSAD-T-CSM 

Fort  Bliss,  Texas  79916 

Commandant  1 

U.S.  Army  Command  and 
General  Staff  College 
Attn:  Acquisition,  Library  Div 
Fort  Leavenworth,  Kansas  66027 

Dr.  Hans  K. Ziegler  (AMSEL-TL-D)  1 
Army  Member,  TAC/JSEP 

U.S.  Army  Electronics  Command  (DRSEL-TL-D) 
Fort  Monmouth,  New  Jersey  07703 

Mr.  J.E.  Teti  (AMSEL-TL-DT)  3 

Executive  Secretary,  TAC/JSEP 
U.S.  Army  Electronics  Command  (DRSEI ,-TL-DT) 
Fort  Monmouth,  New  Jersey  07703 

Director  1 

Night  Vision  Laboratory,  ECOM 
ATTN:  DRSEL-NV-D 

Fort  Belvoir,  Virginia  22060 

Commander/Director  1 

Atmospheric  Sciences  Laboratory  (ECOM) 

Attn:  DRSEL-BL/DD 

White  Sansa  Missile  Range 
New  Mexico  88002 

Director  1 

Electronic  Warfare  Lab.  , ECOM 
Attn:  DRSEI,- WL- MV 

White  Sands  Missile  Range 
New  Mexico  88002 

Commander  1 

US  Army  Armament  Command 
Attn:  DRSAR-RD 

Rock  Island,  Illinois  61201 


-3- 

JSEP  3/77 


No.  of 
Copies 


No.  of 
Cop i es 


Project  Manager  1 

Ballistic  Mi ssi le  Defense  Program 
Office 

Attn:  DACS-BMP  ( Mr.  A.  Gold) 

1300  Wilson  Blvd. 

Washington,  D.C.  22209 

Director,  Division  of  Neuropsychiatry 
Walter  Reed  Army  Institute  1 

of  Research 

Washington,  D.C.  20012 


NI.-H 

Dr. 

F. 

Schwering 

1 

TL-E 

Dr. 

S. 

Kronenberg 

1 

TL-E 

Dr. 

J. 

Kohn 

1 

TL- 1 

Dr. 

C. 

Thornton 

1 

NL-B 

Dr. 

S. 

A mors os 

1 

Col.  Robt.  W.  Noee 
Senior  Standardization  Rep.  1 

U. S.  Army  Standardization 
Group,  Canada 

Canadian  Force  Headquarters 
Ottawa,  Ontario,  Canada  KIA  0K2 


Commander,  USASATCOM  1 

Fort  Monmouth,  New  Jersey  07703 

Commander,  U. S.  Army  1 

Communications  Command 
Attn:  Director,  Advanced  Concepts 
Office 

Fort  Huachuca,  Arizona  85613 

Project  Manager,  ARTADS  1 

EAI  Building 

West  Long  Branch,  N. J.  07764 


Commander 

CCOPS-PD 

Fort  Huachuca,  Arizona  85613 
Attn:  H.A.  Lasitter 

Department  of  the  Navy 

Dr.  Sam  Koslov  1 

ASN  (R&D) 

Room  4E741 
The  Pentagon 
Washington,  D.C.  20350 


U.  S.  Army  White  Sands  Missile  Range 
STEW'S- 1 D-R  1 

Attn:  Commander 
White  Sands  Missle  Range 
New  Mexico  88002 

Mr.  William  f.  Kawai 

U.  S.  Army  R&D  Group  (Far  East)  1 

APO,  San  Francisco,  Ca.  96343 

Director,  TRI-TAC  1 

Attn:  TT-AD  (Mrs.  Briller) 

Fort  Monmouth,  N.  J.  07703 


Office  of  Naval  Research  1 

800  N.  Ouincy  Street 
Arlington,  Virginia  22217 
Attn:  Codes  100 

102 
201 
220 
221 
401 

420 

421 

427  (All  Hands) 

432 

437 


Commander 

U.  S.  Army  Electronics  Command 

Fort  Monmouth,  N.  J.  07703 

Attn  AMSEL-RD-O  (Dr.  W. S. McAfee)  1 


CT-L 

(Dr. 

G.  Buser) 

1 

NL-0 

(Dr. 

11.  S.  Bennett) 

1 

NL-T 

(Mr. 

R. Kulinyi  ) 

1 

TL-B 

1 

VL-D 

1 

W1.-D 

1 

TL-MM 

(Mr. 

Li  pet z) 

1 

( c on  t 

*d) 

Naval  Research  Laboratory 
4555  Overlook  Aven.  SW 
Washington,  D.C.  20375 


Codes  4000  - 

Dr. 

A 

Berman 

4105  - 

Dr. 

S. 

Teitler 

4207  - 

Dr. 

J. 

McCaffrey 

5000  - 

Dr. 

H. 

North 

5200  - 

Mr. 

A. 

Brodz insky 

5203  - 

Dr. 

L. 

Young 

5210  “ 

Dr. 

J. 

Da  vey 

JSEP 


3/77 


No.  of 

No.  of 

Copi es 

Copies 

Naval  Research  Laboratory 

1 

R.  N.  Keeler 

1 

4555  Overlook  Ave.  SW 

NAVMAT  - Code  03T 

Washington,  D.C.  20375 

CP  # 5 

Attn:  Codes  5220  - Mr.  H. 

Lessof f 

2211  Jefferson  Davis  Hwy. 

5230  - Dr.  R. 

Green 

Arlington,  Virginia  20360 

5250  - Cf.  L. 

Whicker 

5260  - Dr.  D. 

Bar  be 

Mel  Nunn 

1 

5270  - Dr.  B. 

McCombe 

NVMAT  0343 

5300  - Dr.  M. 

Skolnik 

CP#  5,  Room  1044 

5403  - Dr.  J. 

Shore 

2211  Jef  erson  Davis  Hwy. 

5464/5410  - Dr.  J.  Davis 

Arlington,  Virginia  20360 

5500  - Dr.  T. 

Jacobs 

5509  - Dr.  T. 

Giallorenzi 

Dr.  F. I . Tanczos 

1 

5510  - Dr.  W. 

Faust 

NAVAIR-03B 

6400  - Dr.  C. 

Klick 

JP#  1,  Room  412 

7701  - Mr.  J. 

Brown 

1411  Jefferson  Davis  Hwy 
Arlington,  Virginia  20360 

Director 

Office  of  Naval  Research 

1 

Dr.  H. J.  Mueller 

1 

495  Summer  Street 

Naval  Air  Systems  Command 

Boston,  Mass.  02210 

Code  310 
JP  # 1 

Director 

1 

1411  Jefferson  Davis  Hwy. 

Office  of  Naval  Research 

Arlington,  Virginia  20360 

New  York  Area  Office 
715  Broadway  5th  Floor 

Mr.  N.  Butler 

1 

New  York,  New  York  10003 

Naval  Electronics  Systems 
Code  304 

Command 

Director  of  Naval  Research 

Branch 

Office 

NC  # 1 

536  South  Clark  Street 

2511  Jefferson  Davis  Hwy. 

Chicago,  Illinois  60605 

1 

Arlington,  Virginia  20360 

Director  of  Navnl  Research 

Branch 

Office 

Mr.  L.  W.  Sumney 

1 

1030  East  Green  Street 

1 

Naval  Electronics  Systems 

Command 

Pasadena,  Calif.  91101 

NC  # 1 

2511  Jefferson  Davis  Hwy. 

Office  of  Naval  Research 

1 

Arlington,  Virginia  20360 

San  Francisco  Area  Office 
760  Market  St.  Room  447 

J.H.  Huth 

1 

San  Francisco,  Calif.  94102 

Harris  B.  Stone  1 

Office  of  Research,  Development,  Test 
Evaluation  NOP-987 
The  Pentagon.  Room  5D760 
Washington,  D.C.  20350 


Dr.  A. L. Slaf kosky 
Code  RD-1 

Headquarters  Marine  Corps 
Washington,  D.C.  20380 


1 


NAVSEA  - Code  03C 
NC  # 3,  Room  11E08 
2531  Jefferson  Davis  Hwy. 
Arlington,  Virginia  20362 

Capt.  R. B.  Meeks 

Naval  Sea  Systems  Command 

NC  #3 

2531  Jefferson  Davis  Hwy. 
Arlington,  Virginia  20362 

-5- 


JSEP  3/77 


No.  of 
Copies 


No.  of 
Cop i es 


Naval  Surface  Weapons  Center  1 

White  Oak 

Silver  Spring,  Maryland  20910 
Attn:  Codes  WR  - 04  - W.  Scanlon 

WR  - 30  - Dr.  J.  Dixon 
WR  - 303  - Dr.  R.  Allgaier 
WR  - 34  - H. R.  Riedl 
WR  - 43  - P.  Wessel 

Naval  Surface  Weapons  Center  1 

Dahlgren,  Virginia  22448 


Attn: 

Codes  DF  - 

J. 

Mills 

DF  - 

14 

- K.  Ferris 

DF  - 

36 

- S.  Leong 

Naval  Air  Development  Center  1 

Johnsvi lie 

Uarminister,  Penna  18974 
Attn:  Codes  01  - Dr.  R.  Lobb 

202  - T.  Shoppie 
20212  - S.  Campagna 
2022  - G.  Fer 

Dr.  Gernot  M.  R.  Winkler  1 

Director,  Time  Service 
U. S.  Naval  Observatory 
Mass.  Ave.  at  34th  St.,  N.  W. 
Washington,  D.C.  20390 

Officer  in  Charge  1 

Carderock  Laboratory 
David  Taylor  Naval  Ship  Research 
and  Development  Center 
Bethesda,  Maryland  20034 

Officer  in  Charge  1 

Annapolis  Laboratory 

Naval  Ship  Research  b Development 

Center 

Annapolis,  Maryland  21402 

Dr.  G.  Gould  Technical  Director  1 
NavalCoastal  Systems  Laboratory 
Panama  City,  Florida  32401 

M.  J.  Wynn  1 

Code  790 

Naval  Coastal  Systems  Laboratory 
Panama  City,  Florida  32401 


Dr.  W.  A.  VonWinkle  1 

Associate  Technical  Director  for 
Technology 

Naval  Underwater  Systems  Center 
New  London,  Connecticut  06320 

Officer  in  Charge  1 

Naval  Underwater  Systems  Center 
Newport,  Rhode  Island  02840 

Dr.  H.  L.  Blood  1 

Technical  Director 
Naval  Undersea  Center 
San  Diego,  Calif.  95152 

Dr.  Robert  R.  Fossur  1 

Dean  of  Research 
Naval  Postgraduate  School 
Monterey,  Calif.  93940 

Naval  Electronics  Laboratory  Center 
271  Catalina  Blvd.  1 

San  Diego,  Calif.  92152 
Attn:  Codes  0220  - H.T.  Mortimer 

2000  - P.C.  Fletcher 
2020  - V. E.  Hildebrand 
2100  - C.A.  Nelson 
2200  - J.  Ritcher 
2300  - C.W'.  Erickson 
2400  - F. M.  Tripak 
2500  - W.E.  Richards 
3000 

3400  - R.  Coburn 
4000  - C.E.  Pierson 
4600  - I.  Lagnado 
5000  - A. E.  Beutel 
5200  - R.R.  Eyres 
5300  - P.H.  Johnson 
5600  - W.  J.  De  jka 

Naval  Weapons  Center  1 

China  Lake,  Calif.  93555 
Attn:  Codes  60  - Royce 

601  - F. C.Essig 

6013  - V. L.  Rehn 

6014  - D. J.  White 

6018  - J. M Bennett 

6019  - N.  Bottka 
(cont ’d) 


-6- 


JSEP  3/77 


No  of 
Cop i es 


No.  of 
Cop i es 


Naval  Weapons  Center 
China  I,ake,  Calif.  93555 
Attn:  Codes  605  - W. S.  Me Ewan  1 

5515  - M.H.  Ritchie 
3945  - D.G.  McCauley 
5525  - Webster 


35  - 

D.  J. 

Russel 1 

55  - 

B.  W. 

Hayes 

3544  - 

H.  W. 

Swi nf  ord 

3815  - 

R.S. 

Hughes 

D. E.  Kirk  1 

Professor  & Chairman,  Electronic 
Eng i neer i ng 
Sp-304 

Naval  Postgraduate  School 
Monterey,  Calif.  93940 

Professor  Sydney  P.  Parker  1 

Electrical  Engineering  Sp-62 
Naval  Postgraduate  School 
Monterey,  Calif.  93940 


Robert  E.  Frischell  1 

Johns  Hopkins  University 
Applied  Physics  Laboratory 
Laurel,  Maryland  20810 

Mr.  G.  H.  Gleissmer  1 

Code  18 

David  Taylor  Naval  Ship  R&D  Center 
Bethesda,  Maryland  20084 

Commander  1 

Pacific  Missile  Test  Center 
Code  4253-3 

Point  Mugu,  Calif.  93042 

Richard  Holden  1 

DE  - 34 

Naval  Surface  Weapons  Center 
Dahlgren  laboratory 
Dahlgren,  Virginia  22448 

Other  Government  Agencies 


Dr.  Roy  F.  Potter  1 

3868  Talbot  Street 

San  Diego,  Calif.  92106 

Mr.  J.C.  French  1 

Electronics  Technology  Division 
National  Bureau  of  Standards 
Washington,  D.C.  20234 

John  L.  Allen  1 

Deputy  Director  (Research  & Advanced 
Technology ) 

ODDR&E 

The  Pentagon,  Room  3E114 
Washington,  D.C.  20301 

Leonard  R.  W’eisberg  1 

Assistant  Director  (Electronics 
& Physical  Sciences) 

ODDR&E 

The  Pentagon 
Washington,  D.C,  20301 


Mr.  F.C.  Schwenk,  RD- V 1 

National  Aeronautics  and 
Space  Administration 
Washington,  D.C.  20546 

Los  Alamos  Scientific  Lab  1 

Attn:  Reports  Library 

P.0.  Box  1663 

Los  Alamos,  New  Mexico  87544 

M.  Zane  Thornton  1 

Deputy  Director, 

Institute  for  Computer 
Sciences  & Technology 
National  Bureau  of  Standards 
Washington,  D.C.  20550 

Director,  Office  of  Postal  1 

Technology  (R&D) 

U. S.  Postal  Service 
11711  Parklawn  Drive 
Rockville,  Maryland  20852 


George  Gamota 

Staff  Specialist  for  Research 
ODDR&F. 

The  Pentagon,  Room  3D1079 
Washington,  D.C.  20301 


NASA  Lewis  Research  Center 
Attn:  Library 

21000  Brookpark  Road 
Cleveland,  Ohio  44135 


-7- 


JSEP  3/77 


1 


I 


No.  of 
Copies 


No.  of 
Copies 


Library  - R51  1 

Bureau  of  Standards 
Ac qui  si t i on 

Boulder,  Colorado  80302 

MIT  Lincoln  Laboratory  1 

Attn:  Library  A-032 

P.0.  Box  73 

I,exington,  Mass.  02173 

Dr,  Jay  Harris  1 

Program  Director,  Devices  and 
Waves  Program 

National  Science  Foundation 
1800  G.  Street 
Washington,  D.C.  20550 

Dr.  Howard  W.  Etzel,  Deputy  Director 
Division  of  Materials  Research  1 

National  Science  Foundation 
1800  G.  Street 
Washington,  D.C.  20550 

Dr.  Dean  Mitchell,  Program  Director 
Solid-State  Physics  1 

Div.  of  Materials  Research 
National  Science 

Non-Government  Agencies 

Director  1 

Research  Lab.  of  electronics 
Massachusetts  Inst,  of  Tech. 
Cambridge,  Mass.  02139 

Director  1 

Microwave  Research  Institute 
Polytechnic  Inst,  of  New  York 
Long  Island  Graduate  Center 
Route  110 

Farmlngdale,  New  York  11735 


Director  1 

Columbia  Radiation  Laboratory 

Department  of  Physics 

Columbia  University 

538  West  120th  Street 

New  York,  New  York  10027 

Director  1 

Electronics  Research  Laboratory 
University  of  California 
Berkeley,  Calif.  94720 

Director  1 

Electronics  Sciences  Laboratory 
University  of  Southern  California 
Los  Angeles,  California  90007 

Director  1 

Electronics  Research  Center 
The  University  of  Texas  at  Austin 
Engineering-Science  Bldg.  112 
Austin,  Texas  87812 

Director  of  Laboratories  1 

Division  of  Engineering  and 

Applied  Physics  - Tech.  Reports 

Collection 

Harvard  University 

Pierce  Hall 

Cambridge,  Massachusetts  02138 


Assistant  Director 
Microwave  Research  Institute 
Polytechnic  Inst,  of  New  York 
333  Jay  Street 
Brooklyn,  New  York  11201 


-8- 


JSEP  3/77 


